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1. Introduction 
We describe singularity formation for the semi-linear wave equation 

(1.1) Ou- ^u{l-u^) = 0, n^dtt-A 

in M2+^. This equation arises as follows: consider Yang-Mills fields in {d + 1)- 
dimensional Minkowski spacetime. The gauge potential Aa is a one-form with 
values in the Lie algebra g of a compact Lie group G. In terms of the curvature 
Faf3 = daAfs — dpAa + [^a, Ap] the Yang Mills equations take the form 

where [•, •] is the Lie bracket on G. We take G = SO{d) with g being the skew- 
symmetric dxd matrices. In particular Aa = Assuming the spherically 
symmetric ansatz (see |14j and [8] for analogous considerations in the context of 
the Yang-Mills heat flow) 

the Yang Mills equations reduce to the semi-linear wave equation 

This equation is invariant under the scaling u{r,t) )■ u{r/X,t/X). With respect to 
this scaling the energy 



E = 



•dr 



is invariant iff d = 4 which is the case we consider in this paper. Equation p.ip 
admits the stationary solution 

1 ^2 

Qir) 



1 + 7-2 ' 

called instanton. In 3 -I- 1 dimensions, the Yang-Mills equations are subcritical 
relative to the energy. Eardley and Moncrief [6] , [7] showed that in that case there 
are global smooth solutions. See also Klainerman, Machedon [9] who lowered the 
regularity assumptions on the data. In the energy critical case of 4 -I- 1 dimensions, 
local well-posedness in with s > 1 was shown by Klainerman, Tataru [lOj . 
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However, it was conjectured that global wellposedness fails and that singularities 
should form, see Bizon, Tabor [1] and Bizon, Ovchinnikov, Sigal [2\ for numerical 
and heuristic arguments to that effect. However, such a phenomenon had not been 
observed rigorously. In this paper we show how to construct a solution to the wave 
equation Ijl.ip as a perturbation of a time-dependent profile 

with X{t) a logarithmic correction to the self-similar ansatz 

\{t) = t-\-\ogtf, /3 > 

In other words, we prove that in general the energy critical Yang-Mills equations 
develop singularities in finite time. 

As in our earlier work [11 for energy critical wave maps, and [12j for the energy 
critical semi-linear wave equation in R"^, the blow up rate is prescribed. Since a 
continuum of rates is admissible, the blow up solutions which we construct can 
of course not be stable. In contrast to the rates X{t) = t~^~'^ which appeared in 
[TT] and [T^], in the case of Yang-Mills we only make logarithmic corrections to 
the self-similar rate. This has to do with the fact that the linearized Yang Mills 
operator has a zero energy eigenvalue in 4-1- 1 dimensions, whereas for wavemaps as 
well as the three-dimensional semi-linear focusing wave equation, it exhibits a zero 
energy resonance. This difference is very important as in the case of an eigenvalue 
an orthogonality condition appears which is not present in the zero energy resonant 
case. It is this fact which required major changes to our scheme, especially to the 
"renormalization part" in which we construct approximate solutions. In addition, 
in contrast to our earlier work on wave maps llj, the approximate solutions here 
are much rougher, and indeed asymptotically only lie in H^, the threshold for local 
well-posedness of the critical Yang-Mills equation. The reason for this is the much 
more singular nature of the ODE's arising in the renormalization step, due to the 
different blow up rate. 

Theorem 1.1. Let 

\{t) :=i-i(-logt)^ 

For each /3 > § there exists a spherically symmetric solution u to (II. 1[) inside the 
cone {r < t,t < to} which has the form 

u{x, t) ^ Q{r\{t)) + v{x, t) 

where the function v has the size and regularity, with S :— tdt + rdr, 

\\Wv\\l2 + \\WSv\\l2 + \\WS\\\l2 < I logir^ 

as well as the pointwise decay 

\v{t,x)\<\logt\-' 

We emphasize that our solutions are just barely better than H^, in contrast to 
our earlier work on wavemaps. While local wellposedness is not known for the 
general Yang-Mills problem, it is known in the equivariant case, see [15]. This is 
important for our purposes, as it shows that our solutions belong to a class for which 
a local wellposedness theory is available. In addition, the vector field S is required 
to control the strong singularity in the nonlinearity r^'^u'^ at r = 0; this is in the 
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spirit of the method of invariant veetor fields in nonhnear wave equations which 
allows for improved decay away from the characteristic light-cone — t}. More 
precisely, one can use elliptic estimates close to r = to control the aforementioned 
singularity. 

2. The proof of the main theorem 

This section contains the proof of Theorem 11.11 The first step is construct an 
arbitrarily good approximate solution to the wave equation (|l.ip as a perturbation 
of a time-dependent profile uq — Q{R). The result is as follows: 

Theorem 2.1. For each integer N there exists a spherically symmetric approximate 
solution UN to (jl.ip inside the cone {r < t,t < to} which has the form 

i?4 



UN{r,t) = Q{R)+vio{r,t)+VNir,t), t;io(r,i) (tXit))- 



4(l+i?2)2 

with X{t) — t^^\ log^l^, R = rX{t), and vj\i satisfying the pointwise bounds 

r2 i?2 



|wjv(r,t)| + \SvN{r,t)\ + \S\Nir,t)\ 



< 



t^\l0gt\ (U(0)2|logt| 

and so that the corresponding error 

2 

satisfies 

|ejv(r,t)| + \SeN{r,t)\ + \S^eN{r,t)\ < t~^\\ogt\-^ 

The proof of the above theorem is carried out in Section [3l The description of 
the approximate solutions upf obtained there is much more precise than what is 
stated above. In particular, the functions itjv are analytic up to the cone t = r, and 
the nature of the singularity at the cone is clearly explained. 

Given the approximate solutions constructed above, we look for a solution u 
to (fLTjl of the form 

u{t,r) = UN{t,r) + e{t,r), 
where e is to be determined via Banach iteration. The equation for e is 

/ 1 2 \ 

{ -dt +d^ + -dr + ^(1 - 3uNit,rf - 3e{t,r)uNit,r) - e^{t,r))]e{t,r) = cn 

We divide this equation into a linear part and a nonlinear perturbative term. Based 
on past experience one would expect that in the main linear part M^r is simply 
replaced by Q{X(t)r). However, in this case that is not enough. Instead, as it turns 
out, some of the effects of the first correction term wio also need to be taken into 
account. Hence the above equation is rewritten in the form 

(2.1) f - 92 + dl + -dr + ^{l- 3Q(A(t)r)2 - &Q{X{t)r)vw)]e - ejv + AA(e) 
where 

M{e) = ^ (3£« - QiX{t)rf - 2Q{X{t)r)vio) + Ss^un + 
We first consider the linear problem 

(2.2) (-d^ + dl + -dr + ^(1 - iQ{X{t)rf - QQ{X{t)r)vw))e = f 
\ r r J 



4 



J. KRIEGER, W. SCHLAG, AND D. TATARU 



where the principal spatial part is given by the selfadjoint operator 

Lt = -dl - ^dr - ^(1 - 3Q{X{t)r)^) 

This is time dependent, but is obtained by rescaling from the operator 

L = -^(l-3Q(r)2) 

We remark that, as proved in the next section, L is a nonnegative operator. 

A difficulty that we face in solving (|2.ip iteratively is in handling the singularity 
at in the e'^ term in AA(e) . Energy estimates on e do not suffice, so we introduce 
the scaling vector field 

S tdt + rdr 

and we seek to simultaneously bound e, Se and S'^e in a norm that is a scale 
adapted version of the norm, 

\\4Hh — sup llogtl^-f^-Wm e\\L^rdr) + \\dt£\\L^{rdr)+ H't)\^0gt\-l^\\e\\L2(rdr)) 
Q<t<ta 

For /, on the other hand, we just use uniform bounds, 

||/|U.^:= sup \-\t)\\0gtf\\f{t)\\L2^rdr) 
0<t<ta 

The main result of the linear theory is the following theorem. It is proved in 
Section [6l 

Theorem 2.2. There exists a linear operator $, so that for each f the function 
e = $/ solves (12. 2|) . and for all large enough Nq 3> Ni 3> N2 it satisfies the bounds 

The implicit constants here depend only on (3. 

We note that $ is in effect the forward solution operator for the equation (|2.2p . 
In this theorem / is not required to be supported inside the cone {r < t}. However, 
if that is the case the $/ is also supported inside the cone due to the finite speed 
of propagation. 

In order to prove the above theorem it is convenient to pass to different coordi- 
nates in which the Schrodinger operator is no longer time dependent. Specifically, 
introduce new coordinates (r, R) given by 

i-to 

X{s)ds, R = A(t)r 

Jt 

Then, denoting 

i{T,R) ■.^R^e{t,r), f{T,R) -.^ R^ f [t , r) 
where A is now understood as a function of r, the equation (|2.2p becomes 

-{dr + ^Rdnf + \{^f + \dr{^)\ e- CI- ^^Q{R)v,,I= /, 



(2.3) 




(2.4) 




(2.5) 
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where 

_^ 24 

dm ^ 14i?2 " (l + i?2)2 

The spectral properties of the operator C are studied in Section [H These are 
essential in the proof of Theorem 12.21 in Section [6l 

We next continue the proof of our main result with the perturbative argument 
for the equation (|2.ip . By the construction in Section [3] we know that for arbi- 
trarily large A^o 3> Ni 3> N2 we can find an approximate solution un so that the 
corresponding error cn satisfies 

l|eAr||y := \\eN\\L%^_^ + ll^eArH^^^^ + \\S\n\\lI^ « ^ 

where the smallness is gained by taking to small enough. It is important to note 
that, even though cn has limited regularity, the roughness is relative to the self- 
similar variable a := ^ which satisfies Sa = 0. For this reason S^cn does not lose 
any regularity. We iteratively construct the sequence {{ej, fj)}j>o by 

fo:^eN, Ej-.^^fj, fj+i := CN +Af{ej) 

and show that it converges to a solution e of (|2.1[) in the norm 

+ \\Se\\„.^ + \\Sh\\H^^ 

Nq Ni N2 

By Theorem 12 . 2 1 we know that $ is a bounded operator with smah norm, 
(2.6) ll*/llx<iVo"'ll/lli' 

The proof is concluded if we show that the nonlinear term satisfies a similar bound: 

Lemma 2.3. The map f J\f{^f) is locally Lipschitz from Y to Y , with Lipschitz 
constant of size 0(iV^ ). 

Proof. We denote e = <&/, and successively consider the linear and the nonlinear 
terms in A/'(e). 

A. The linear term has the form 

Ui{e) = ge, g = ^« - Q{X{t)r)^ - 2Q(A(<)r)z;io) 
By construction we have 

l5l + 1^51 + < - X{tf\\ogtr'^~' 

which directly leads to 

\me)\\Y<\\e\\x. 

where only the components of the H]^, norms are being used (as part of the 
space X). The desired conclusion now follows from (|2.6p . 

B. The nonUnear term has the form 

A/'2(e) = ^{3uNe^ + e^) 

The coefficient un satisfies 

\un\ + \Sun\ + \S\n\<1 
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SO we can neglect it. The main difficulty here arises from the singular factor on 
the right-hand side. To address that we will establish several bounds. The first is 
a pointwise bound, 

(2.7) H<|logir+'''-^||HlHj, 
which applied to e, Se and yields 

(2.8) <_L|iogi|i+2/3-^^||/|l^, fc = 0,l,2 

JVfe 



The second is a weighted bound, namely 

(2.9) \\r-\{t)\\L. <\{t)\\ogt\-^^+^\\f\\y 
Interpolating between the k — 2 case of (|2.8p and (|2.9p we also obtain 

(2.10) \\r-'Se{t)U^ < ^A(t)^ | logr'^^+'+''ll/||y 



By ((2?8)) with k = and ((2^ we obtain 
(2.11) 



L2 



<-lA(t)|logtp+2^-^o-^.(|,;|,^ + |,;|,^3) 



(2.12) ||5AAi(e)(t)||^, <_A(<)|logtp+2^-^-^^(||/||2, + 11/11^3) 



Using also (12. Sp with fc = 1 we similarly obtain 

1 

^NO ^ Nl' 

Finally, due to (|2.8p with k = 2 and (|2.10p we also have 

(2.13) ||52AAi(e)(t)||^, K^xmogtl^+^^-^^^iWffy + WfWrs) 

Together, the bounds ()2.1ip - (|2.13p suffice to obtain the conclusion of the lemma 
provided that A''2 is large enough. 

It remains to prove the bounds (|2.7p and (12. 9p . For the operator L we have the 
straightforward elliptic bound 

||Vu;|U2 + \\r-^w\\L2 < \\Liw\\L2 + \\w\\l2 
By rescaling this gives 

.2 \\^M\l2 + \\r-'wh2 < \\LfwU2 + A(t)||u;|U. 

Then (|2.7p follows from the point-wise bound for spherically symmetric functions 
in R2 

I|w||l- < ||Vm||l2 + \\r-^u\\L2 
Next we turn our attention to the bound (|2.9p . Due to (|2.8p (fc — 0) it suffices 
to consider the region r < t/2. In this region we use the scaling vector field 
S — tdt + rdr to derive a stronger equation for e. From 

tdt — S ~ rdr 

one infers that 

t^d^e + tdte = -S^e + r^d^e + 2tdtSe 
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and further 



dr 



1 



-dr - — e = t-^{-S'^e + 2tdtSe - tdte) - (□ + — )e 



Due to 



r 



we can estimate the last term by 
4 



□ 



< 



\f\ + y\e\ 



This leads to the bound 



< i-2[||^2^||^, + \\{tdtS)s\\L^ + \\{tdt)£\\LA + 



< 



+ A'lkl|L^ + ll/l|L^ 



X{t)\logt 



,f|l-JV2| 



l/ll 



L2 



Taking also into account (|2.14p applied to e with N = Nq, the estimate (|2.9|) would 
follow from the fixed time bound 



(2.15) \\r-^e\\L. <\\t-'^eU2 + \\t-'r'h\\L. 



- r 



■dr 



1 



,9 ^'rr dr ^ I £ 

r 



L2 



This rescales to t = 1, in which case it is a standard local elliptic estimate near 
r = 0. □ 



3. The renormalization step 

In this section, roughly following [KSTl], we show how to construct an arbitrarily 
good approximate solution to the wave equation p.ip as a perturbation of a time- 
dependent profile 

(3.1) uo = Q(i?), R = rX{t), $(i?) = _J^ 

with X{t) a logarithmic correction to the self-similar ansatz 

\{t) = t-\-\ogtf,f3>l 

In fact, for ease of notation we will take f3 £ Z; the general case is only a minor 
modification of the integral one and we leave it to the reader. This ansatz is quite 
natural in light of a necessary orthogonality condition which makes its appearance 
in the ensuing considerations. We note, however, that by contrast to [KSTl], the 
approximate solutions here are much rougher, and indeed asymptotically only lie 
in H^, the threshold for local well-posedness of the critical Yang-Mills equation. 
The reason for this is the much more singular nature of the ODE's arising in the 
renormalization step, due to the different blow up rate. 

The following is the main theorem of the first half of the paper. Throughout this 
section, we will work on the light-cone {r < t} (in particular, all functions in this 
section will be defined only on r < i). 

Theorem 3.1. Let k G N. There exists an approximate solution U2k-i G for 
(|l.ip of the asymptotic form (as R oo) 

u,,Mr,t) ^ QiXit)r) + ^^^^-^ + O ((^xy^j 
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SO that the corresponding error has size 

e2k-i = O 



t^{tX) 



2k 



Here the O(-) terms are uniform in < r < t and < t < to where to is a fixed 
small constant; they are also stable with respect to the application of powers of the 
scaling operator S. We also have M2fc-i(-,i) € C°°{[0,t)), and further U2k-i G ■ 
The only singularity arises on the light cone r — t. 

Proof. We iteratively construct a sequence Uk of better approximate solutions by 
adding corrections Vk, 

Uk=Vk+ Uk-l 

The error at step k is 

-drjUk - ;^ 



Ck = (df - dl - -dr)Uk - -f{Uk), f{u) - U{\ - u") 



To construct the increments Vk we first make a heuristic analysis. If u were an 
exact solution, then the difference 

£ = w - Uk-i 

would solve the equation 

{-d^^+dl + Ur)e+^f'{uk-i)e 
2 

(3.2) = gfc-i - ^(/(^ife-i + e) - f'{uk-i)e - /(ufe_i)) 

2 

= Ck-i + -^{ie^uk-i+e^) 

In a first approximation we linearize this equation around e = and substitute 
Uk-i by Q{R). Then wc obtain the linear approximate equation 

(3.3) (-d^ + dl + -dr + 4(1 - 3g')) e « Ck-i 

For r <C t we expect the time derivative to play a lesser role so we neglect it and 
we are left with an elliptic equation with respect to the variable r, 

(3.4) (^ol + ldr + ^{l-3Q^)^e^ek-i, r^t 
For r K t we rewrite (13.31) in the form 



1/) 4 

Here the time and spatial derivatives have the same strength. However, we can 
identify another principal variable, namely a = r/t and think of e as a function 
of {t,a). As it turns out, neglecting a "higher order" part of e^^i which can be 
directly included in Ck, we are able to use scaling and the exact structure of the 
principal part of ek-i to reduce the above equation to a Sturm-Liouville problem 
in a which becomes singular at a = 1. 
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The above heuristics lead us to a two step iterative construction of the v/c's. The 
two steps successively improve the error in the two regions r <^ t, respectively r k t. 
To be precise, we define Vk by 

(3.5) (dl + -dr + \(l- 302)") = e°,. 
respectively 

(3.6) {-d^, + dl + -dr - 4) ^'2fe = 

\ r r'^ J 

both equations having zero Cauchy datsQ at r = 0. Here at each stage the error 
term Cfe is split into a principal part and a higher order term (to be made precise 
below) , 

The successive errors are then computed as 

e2k = e2/c-i + N2k{v2k), e2fe+i = e\k + d1v2k+i + A^2fc+i(w2/c+i) 

where 

(3.7) N2k+i{v) = ^{ult, - Q'')v + ^{3v\2k + v^) 
respectively 

(3.8) N2k{v) = ^{ujk-i - 1)« + ^{3v^U2k-i + v^) 

To formalize this scheme we need to introduce suitable function spaces in the 
cone 

Co = {{r, t) : <r <t,0 <t <to} 

for the successive corrections and errors. We first consider the a dependence. For 
the corrections Vk we use the following general concept 

Definition 3.2. Let k > 0. Then Qk is the algebra of continuous functions 
q{a, a, ai) 

g : (0, 1] X M X R -> M 

with the following properties: 

(i) q is smooth in a £ (0, 1), and meromorphic and even around a — 0. Further, 
the restriction to the diagonal 

q{a, b) := q{a, b, b) 

extends analytically to a = and has an even expansion there. 

(ii) q has the form 

j<0,i<b|/2 

q{a,a,ai)= ^ (7y(a,loga,logai)a'aj 

i+j<-k 

with qij polynomial in log a, log ai. The sum only has finitely many terms. 
(Hi) Near a = 1 we have a representation of the form 

q = qoia,a) + (1 - a^)^qi{a,a,ai) 
with coefficients qo, qi analytic in a around a = 1. 



^The coefRcients are singular at r = 0, therefore this has to be given a suitable interpretation 
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The order of the pole at a = as it appears in Definition 13.21 part (i), is 
controUed by some absolute constant depending only on k. The same comment 
applies to every pole at a = appearing in this section and will be assumed tacitly 
throughout. For the errors we introduce another functions class: 

Definition 3.3. Qj. is the space of continuous functions q{a, a, ai) 

g : (0, 1) X K X M ^ M 

with the following properties: 

(i) q is smooth in a ^ (0, 1), meromorphic and even around a = 0. The restric- 
tion to the diagonal 

q{a, b, b) 

extends analytically to a = 

(a) q has a representation as in (ii) of the preceding definition 
(Hi) Near a = 1 we have a representation of the form 

q = qo{a, a) + (1 - a'^)iqi{a, a, ai) + (1 - 0^)^^(72(0, ai) 

with coefficients qo, qi, q2 analytic with respect to a around a — 1. Moreover, 52 
has the same representation as q in (ii), but with k replaced by k + 1 and j < —1. 

Next we define the class of functions of R: 

Definition 3.4. S"^ {R'^ {log RY) is the class of analytic functions v : [0, 00) M 
with the following properties: 

(i) V vanishes of order m at R — 0, and R^"^v has an even expansion around 
R = 0. 

(ii) V has a convergent expansion near R = 00, 

v{R)^ c,,R''-^\\ogRy 
o<j<e+i 

Finally, we introduce the auxiliary variables 

5:-|logi|, 61 := |logi| + |logp(a)| 

where p is a real even polynomial with the following properties: 

j5(l) = 0, p'(l) = -l, p{a) = l + 0{a^') as a^O 

where M is a very large number (depending on the number k of steps in our 
iteration), and p has no zeroes in (0, 1). We can now define the main function class 
for our construction. 

Definition 3.5. a) 5'™(i?'^(log i?)^, Q„) is the class of analytic functions 

V : [Q,oo) X [0,1] X ~> M 

(i) V is analytic as a function of R, 

V : [0, 00) — > Qn 

(ii) V vanishes of order m at R — 

(Hi) V has a convergent expansion at R — 00, 

v{R,;bM)= C,,{;b,b,)R''-^\\0gRy 
0<j<£+i 

with coefficients Cij G Qn- 
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b) I S™' {R'' (iog RY , Qn) is the class of analytic functions w on the cone Co which 
can be represented as 

w{r, t) = viR, a, b, bi), v G 5™(i?'=(log Q„) 

We note that the representation of functions on the cone as in part (b) is in 
general not unique since R, a, b are dependent variables. Later we shall exploit this 
fact and switch from one representation to another as needed. We start our con- 
struction with some explicit computations which allow us to establish the regularity 
of the first few terms in the iteration, namely 

(3.10) t'e, e {Xty' (iSHl) + jr^IS*{R') + -1—IS\R') 

(3.11) «2 ea^/^(l,Qi) 

After these few steps we reach the general pattern, and prove by induction that the 
successive corrections Vk and the corresponding error terms Ck can be chosen with 
the following properties: 

(3.12) V2k-i e I S\R^ [log rY-\ Q2f3k) 

(3.13) t^e2k-i e I S^R^ {log R)''-\ Q'^p^) 

(3.14) V2k G aHS{{\ogRf-\ Q^pik-i)) 

(3.15) t'e2k e a'lSiilogRf-\ Q^^,) + IS\{logRf-\ Q2pk) 

The properties (|3.9p - p.l5p suffice in order to reach the conclusion of the theorem. 
We note that is easy to verify that all the above classes of functions arc left invariant 
by the scaling operator S. 

The proof of (f3Jl) - ([3T5l) roughly follows that in [WM], [SL]. There is, however, 
an important difference near the light cone: for the critical Wave Maps problem as 
well as the critical focussing semilinear equation, the singularity at the boundary 
of the light cone is well modeled by the expression (1 — a)2+'^, which comes from 
the choice of blow up rate t~'^~^ . For Yang-Mills, due to the much faster blow up 
speed, we need to essentially use the much more singular expression 

(3.16) (^-"^^ 



|logi|-f |logp(a)| 

where p{a) is a polynomial so that p(l) = 0. This renders the algebra significantly 
more delicate. We remark that (j3.16p appears canonically in this section. On one 
hand, (1 — 0)2 is part of a fundamental system of that ODE which ()3.6p reduces 
to in the self-similar variable a = This is exactly what one would obtain by 
neglecting all but the selfsimilar components of the wave operator. However, unlike 
in [WM] , [SL] , here we encounter a nontrivial non-selfsimilar effect which forces the 
exact denominator in (|3.16p . In particular this saves the day by insuring that p.l6p 
belongs to i?^(0, 1) which of course is a minimal requirement here. 
To commence the construction of the Vk, we recall that 
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where $ is the zero eigenfunction, L<i> = with 
By (|3.ip we have 



Step 1: Begin with cq as above and choose vi so that (I3.12p for k = 1 holds. 
Further the error ei thereby generated is of the form (j3.13p for k ~ I. 
Here, we simply put := cq. Reformulate the equation for vi as follows: 

{tXfZVRv.^VRt'eo, L = dl-^ + j^^^ 

Using the above calculation of bq, it is then straightforward to write down an 
absolutely convergent Taylor expansion of vi around R = 0. Since t^eo vanishes of 
second order at 0, it follows that vi vanishes of order four at 0. 

Now we turn to the expansion around R = oo. The leading term in i^eo is 
+ A key fact is that this satisfies the orthogonality condition 

(i?$'(i?) + <I>(i?),$)R2 =0 

It is partly this orthogonality condition which motivates our choice of X{t). As a 
consequence, the solution to Lvio — + does not grow at oo, precisely 

it equals 



Via = -r{tX) 



4' ' (l + i?2)2 

For the remaining terms we do not have such a precise representation since we lack 
the orthogonality condition. We use this fundamental system of solutions for L: 

, i?3 -l-8i?2 + 24i?4logi? + 8i?6 + i?« 

w(-R) = , BONO , ^a\R) 



Their Wronskian is M^((?!)o, 6*0) = 1. Then R-^(f>a{R) and define e(i?) 

R^^OoiR) so that L$ = 0, = 0, respectiveljQ. One thus obtains an integral 
representation for vi using the variation of parameters formula, which gives 

(tXfviiR) = rMR')VRt^ea{R') dR' - [\{R')VRt^eoiR') dR' 

VR Jo vR Ji 

= Q{R)j ^{R')t^eo{R')R' dR' -^{R) J e{R')t^ea{R')R' dR' 

In the end we obtain the representation 

(3.17) VI = «io + vn, vn e (Xt)-^ (^—IS\R^) + .^—IS^{R^ 

V|logi| |logt|2 

which implies (|3.9p . 

Next, we determine the error, which is given by 

t^ei=t^dfvi-^-^{3vlQ + vl) 



Note the appearance of "I'(-R) logi?, as part of 0. 
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After some computations we obtain the relation ()3.10|) . namely 

(3.18) t'e^ e {Xtr' (lS\l) + ^IS*{R') + -1--ISHr' 

Step 2: Recall that V2 is determined by p.6p . which requires specifying e?. This 
will be done iteratively, which means that 

(3.19) e? = X:e?^ 

where J = J{l3) grows with (3 and e'j'"' is specified recursively. To being with, we 
extract the leading order (in terms of growth in R) from ei and set 

with suitable constants ci,C2. Note that then 

ejo :^ei-efe 752(1,22/3) 

which is admissible for 62, see p.lSp . Replacing Q by 1 we now seek to solve the 
linear differential equation 

(3.20) (^-d^ + dl + Ur - ^) = ef 

In the a, b coordinates the above equation is rewritten as 

? 9 

LabV2(a, 0) = Ci — + C2J-^ 

where 

Lab =-{db + adaf - {db + ada) +dl + -da-^ 

Set also the b independent part 

La = il- a^)dl +(-- 2a) da - \ 

For technical reasons, we will only obtain an approximate solution V2 of (|3.20p . We 
then face a dichotomy: either the error Lab W2 — ej" is acceptable for 62 or not; in the 
latter case, we repeat the procedure by including the unacceptable error in and 
solving for a correction to V2- This process (which also needs to take the nonlinear 
component of 62 into account, see (j3.8p ) then leads to the aforementioned iterative 
construction of Cj* and V2- 

2 

We begin by constructing an approximate solution to Lab ^2 = \- The approx- 
imate solution in the following lemma is called W2 rather than V2 since the latter 
will be the sum of various expressions, the first being W2- 

Lemma 3.6. Let e(a) be even analytic and quadratic at a = 0. There is an 
approximate solution W2 for 

LabW2 = b~^e{a) 

which is of the form 

(3.21) W2{a,b)^b-^W^{a)+b^^{l~a)^Wl{a) 
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where W2 , W2 are even analytic in a G (0, 1] with an leading term at so that 
W2 vanishes to fourth order at a = 0. The error has the form 

/a := Lab W2 - 6~^e(a) 

(3.22) = E^{a, 6-1) + (1 - a)iEl{a, b^') 

where E^, E\ are analytic in a € (0,1], linear in b~'^, b^^, b^^, b^^, and vanish 
quadratically at a = 0. 

Proof. We begin with the ansatz 

W^ja) (l-aj^irj(a) 

- + b[ 

where 

LaW°ia) = e(a), La((l - a)^W^{a)) = 
The solvabihty of these equations will bo discussed later in the proof. Then 

KW°{a) K{{l-a)iW2\a)) , e{a) , 

LabW2 = ^ + +/2= — +/2 

where 

/o =i-d^ - ^dbada - dt,)^^^ - (1 - af^Wi{a){d^ + 9(,)[^] 
+ (1 - a)5W^i(a)((a-i - 2a + l)^^ - (1 - afdl)[h 
+ 2(1 - a')da{{l - a)--Wi{a))da[l ] - 2a(l - a)haWi{a)db[^] 

01 Ol 

- 2aWi{a)dadb[^-^^^^] + (1 - af^Wiia)i-da + ((1 - af + (1 - a^))dl)[h 
bi Ol 

The final term here is the same as 

{l-a)iw^{a){-da + 2{l-a)dl)[y^] 

= W^ia)i-da + 2(1 - a)dl)[^^^^] + 2(1 - a)^ Wi(a)9„[^ J 
which implies that the error equals 

/o = i-dl - 2dbada - 96)^^ - (1 - a)iwi{a){dl + db)[y^] 
+ (1 - a)^wUa){{a-' - 2a+ 1)5„ - (1 - a^dDih 

Ol 

(3.23) + (1 - a)i{W^{a) + 2(1 + a)daW^ {a))da[l ^] 
-2a{l-a)idaW^{a)db[^] 

Ol 

(3.24) + W^{a){-2adbda + 2(1 - a)dl - 5„)[-^^ " ' 1 



bi 



RENORMALIZATION AND BLOW UP FOR THE CRITICAL YANG-MILLS PROBLEM. 15 



In the first term we gain at least one power of 6 ^. In the second and fifth terms 
we gain at least one power of bi^. Since 

^^ \a h (1 - a)p'{a) 

(1 - a)dabi = 

p{a) 

which is analytic in [0, 1] it follows that in the third and fourth terms we gain at 
least one power of without losing any power of (1 — a). 

So far we have considered the negligible terms. The key expression is the one in 
the final term, which determines the choice of our ansatz. Here there is a nontrivial 
cancellation which yields an additional 1 — a factor. To begin with, recall that 

{2{l~a)dl~da){l-a)^ =0 

This implies that in (|3.24p at least one derivative has to fall on hi leading to a gain 
of at least one power of bi. However, we need to check that there is no loss in terms 
of powers of (1 — a). This can be seen via the factorization (we first consider dadb 
since the difference from adadt gains a factor of 1 — a) 

{^2di,da + 2(1 - a)dl - da)il - a)^g(a, b) 
(6.25) , , 

= (2(1 - a)-2da - (1 - a)-^){ -db + il- a)da)gia,b) 

provided g{a,b) is smooth. In particular, setting g{a,b) — 
i-2dhda + 2(1 - a)dl - 9,) 

- (2(1 - a)^da - (1 - a)-^){ ~ db + {l- a)da)[h 



- (2(1 - a)ida - (1 - a)-^)6-2(4 ^ (i _ a)da)bi 
(3.26) = (2^^T^5. - 4^^^a.6, - ii^) {d, - (1 - a)0^)b. 
Given our choice of bi 



i^db + (1 - a)da)b, = -1 - 44(1 - a) = 0(1 - a) 

Thus, the (1 — a)-gain in the second factor in p.26p cancels the (1 — a)-loss that 
we incur in the first factor. At the same time we get at least a bi^ factor. In 
conclusion, 



.2 o N(l-a)" ^/(l 



(3.27) {-2adbda + 2(1 - a)d', - Oa) ' ' = O C ,2 ' ) 

where the 0(-)-term here depends linearly on b^^ and b]\^. This establishes the 
desired estimate on the error /2 . 

We now consider the principal part, for which we need to solve 

(3.28) LaW^ = eia), LaHl ~ a)iw^ (a)) = 

In order to analyze these equations, we first discuss fundamental systems of La and 
their respective behaviors at the regular singular points a = and a = 1 of La (we 
can ignore the regular singular point a = — 1 of ia). From 
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we conclude that La[a*^(l + a^0±(a))] ~ where (\)± are even analytic func- 
tions around a = 0. Moreover, a particular solution to La{f) — is given by 
/(a) = — Similarly, for any e(a) as in the statement of the lemma there is 
a particular solution /(a) to La f = e with / even analytic around a — and 
vanishing quadratically at a = 0. Note that / is not unique. However, adding a 
suitable multiple of the a^-homogeneous solution we can achieve that /(a) vanishes 
to fourth order at a = (i.e. /(a) = 0(a^)) and is unique. 
To analyze a fundamental system around a = 1 we write 

La - 2(1 - a)^da{{l - a)^da) - (1 - afdl + a'^l + 2a)(l - a)da - \ 

— '■ Lafi + La.l 

where ia.o '■= 2(1 — a)^da{{^ — a)^da)- Now 

La^oi^ - ^ k{2k - - af-^ 
La{l - af = fc(2fc - 1)(1 - af-^ + 0((1 - af) 
with an analytic 0(-)-term. This implies that Laipo — Latpi = with 
(3.29) ^Ao(a) - l + (l-a)^o(a), i/'i(a) - (1 - a)^(l + (1 - a)V^i(a)) 

where ipa, Vi ^-re analytic around a = 1. In particular, we can solve for in (jX^ 
and W2 is unique up to a constant factor. For future reference we remark that 



4 

La = Plda{p2da) 5- 

pi(a) = -V^ - a^, P2{a) = a\/l - a? 



To solve (jX^ . we first solve for 11^ := 11^2° + (1 ~ a)iw^ and then extract M/j" and 
W2 from it. The logic here is as follows: At a = we want W2 to vanish to fourth 
order. This implies that W must also vanish to the same order since 

5i - 6 = I logp(a)| = - log(l - 0{a'')) = 0{a^') 

with M large. Therefore, as discussed above, W is uniquely determined as a solution 
to 

LaW{a)^a^, -£ < a < 1 

where e > is some small constant. By variation of parameters there exist unique 
constants cq, ci, C2 with the property that 

W^(a) = co?/'o(a) + ciV'i(a) + C2 / [?/'o(a)i/'i(") - V'i(a)V'o(u)](pi(w))~^w^ 



By inspection, the integral on the right-hand side is smooth around a = 1. This 
shows that we need to set 

W^{a) :=ci7/;i(a) 

14^2 («) co'(/'o(a) + C2 / [ip(i{a)%l)i{u) - 7pi{a)Tpo{u)]{pi{u)y^u'^ du 



Observe that at a = we have no guarantee that i ^2 ^''^ smooth; in fact, they 
may exhibit a"^-type behavior. □ 
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(3.30) 



We remind the reader that bi in (I3.2ip cannot be replaced with b since we require 
that W2 G ^^^(0, 1) relative to the a variable. The proof also shows that one cannot 
dispense with the (1 — a)' part of W2 since it is part of the fundamental system 
of La- Another important feature of the previous proof is the cancellation in (j3.25ll . 
For our purposes, g{a,b) = h{bi) whence (I3.25|) becomes 

i-2dkda + 2(1 - a)dl - da){l - a)ih{b,) 

= (2(1 - a)ida - (1 - a)--^)h\bi){- 9b + (1 - a)da)bi 

= (2(1 - a)^h"{bi)dabi + 2(1 - a)^h'(bi)da - /i'(6i)(l - a)-^)0(l - o) 

= 0((1 - a)U"(6i)) + 0((1 - a)U'(6i)) 

In view of (|3.30p . the proof of Lemma [3^ generalizes to right-hand sides such as 
for any fc > 1. 

If we were to now set V2 ■= W2 (from the previous lemma), then the error /2 
from ()3.22p . as well as the remaining C2p" piece from e™, would have to be included 
in 62- However, if /? > 1 this is inadmissible since the error 62 needs to decay at 
least like (fA(t))^^ = b^'^^ . The importance of {tX)~^ lies with scaling; indeed, the 
elliptic equation (|3.5p scales like which equals (iA)^ at its largest. 

These are not the only obstacles we face here: the nonlinear part of 62 (again if 
V2 = ■W2) is 

(3.31) ^{ul~l)w2 + ^{3wlui + wl), 

where ui = Q + vi, see p.Sp . One easily checks that the preceding expression times 
lies in 

(Xty^IS^l, Qi) + iXt)-^IS%R^ Q2). 
The term {Xt)^'^IS'^{l, Qi) can be incorporated into t^e2; however, the term 

iXtr^IS*iR^,Q2) 

is not acceptable for 62 due to the R^ growth. 

We deal with these obstacles by including all unacceptable errors e (with regard 
to 62) in el and solving LabW = e. For example, using the notation of Lemma 13.61 
the second term in p.3ip contributes 

_^ {l-aY^W^{a)W^{a) 

where we replaced ui with 1. The corresponding ansatz for w would then necessarily 
contain the term 

{bbi)-^{l-a)^W{a) 

If dadh (which is part of Lab) hits this term, then we obtain (amongst others) the 
error term 

(l-a)-55-i&r' 

Iterating once more with this error on the right-hand side produces the expression 

(1 - a)^b^'^\ogb 

In order to remove possible singularities at a = (as in the previous proof) one needs 
as many powers of log 61 as of log b. These observations should serve to motivate 
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the following result which will finally allow us to carry out the full iteration for V2 
(as well as for V2k in Step 4 below). We begin with a definition. 

Definition 3.7. Let 2k> m> k> 1. By Tk,m we mean the function class 

m 

Tk,m ■■= {h I h = 6-'=eo(a,log6) + (1 - a)^ ^ e^"(a, log6, log6i)6^-'=6^^' 
+ (1 - a)-5 ^e](a,log6,log6i)6^-'=-i6r^ } 

where for each j the functions eo, e°(a), ej(a) are smooth in a G (0,1), analytic 

around a = I, meromorphic and even around a = 0. Moreover, these functions are 
polynomials in the variables log 6, and logbi, respectively. Further, fk = O(a^) as 

Recall that the order of the pole at a is controlled by a constant depending 
only on k. In what follows, we will tacitly assume that M in the definition of bi 
is sufficiently large depending on k (in fact, the order of the pole at a = in the 
previous definition). Since we are only going to consider finitely many k, this is not 
an issue. Since log6i — log 6 = O(a^) we see that fk{a) = O(a^) is therefore the 
same as 

m m 

eo (a, log 6) + ^ ( 1 - a) ^ e° (a, log 6, log &) + ( 1 - a) ~ ^ ^ fo- ^ (a, log 6, log &) = O (a^ ) 
j=i j=i 

The left-hand side is a polynomial in log b, b~^, so this amounts to the corresponding 
condition for each of its coeflicients. Now for the main iterative lemma. 

Lemma 3.8. The equation 

(3.32) LabV = fk€j^k,m 

admits an approximate solution 

m 

v{a,b) = b-''Vo{a,logb) + (1 - a)5 ^y,(a,log6,log6i)6^-'=6r' 

where Vq, Vj are smooth in a G (0, 1), analytic around a = 1, meromorphic around 
a = 0, and polynomial in the variables log b, log bi . Moreover, v vanishes to fourth 
order at a = and 

Lab V — fk & ^k+l,m + ^k+2,m 

Proof. Let £{0) be the order of the polynomials appearing in the definition of fk 
relative to log 6, and the order relative to log6i with 1 < j < m. We first 
re- write the source term: choose a smooth partition of unity (f)i,2{a), subordinate 
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to the cover (0, 1) — (0, 2e) U (e, 1) for some small £ > 0. Then write 

m 

0i(a)r'=eo(a,log6) + (1 - a)^ ^ e^"(a,log&,log6i)6^-'=6^^' 

711 

+ (1 - a)-^J2e]{a,logb,logb,)V-''-H;'] 

m 

(3.33) = (a) [b-^eoia, log 6) + (1 - a) ^ ^ e°(a, log b, log 

(3.34) + (1 - a)-^ ^e] (a, log 6, log 6)6-^-1] 

(3.35) + (logb- log6i)0i(a)/fc + {b~ bi)Ma)9k+i 

where fk, gk have the same properties as fk- Note that in the expression in brackets 
in p.33p and (|3.34p . all singular powers cancel. For p.35p . expand 

1 I / \| ^ / log \p{a) \ \-i 

01 (a) [log fei - log 6] = (/.i(a)log(l - °^'f°^' ) = -0i(a) ^_^] + error 

Here we may achieve arbitrarily fast decay in time for the error term upon choosing 
N large enough, and hence we can discard its contribution. However, now all the 
terms in 

0i(a)(logb(a)|)-'"/fe, (/.i(a)(logb(a)|)^g,+i, j > 1, 
are smooth up to a = 0, and so are all terms in 

m 

02(a)A - 02 (a) [ 6-^60(0, log 6) + (1 - a)5 ^eO(a,log6,log6i)6^-'=67^' 

m 

+ (1 - a)-5 ej(a, log&, \ogbx)V-^-^b-^\ 
j=i 

These considerations show that we may as well assume that cq, e^, ej are each 
analytic at a = as well as of the form O(a^). With v as in the statement of the 
lemma, we compute 

m 

LaiV= b-''L,Vo{a,\ogb) + Y,V-%^'La{{l - a)^ VS(a, log 6, log6i)) 

rn 

+ Yadb{V^^bl\l - a)"^V,(a,log&,log6i)) + error 

where bi is treated as a parameter, i.e., no derivatives fall on it. Here the last term 
comes from dadb in Lab with the da applied to (1 — a) 2 and db applied to b or log b. 
Assuming that Vj are smooth and that v vanishes of order four at a = one sees 
that the error has the desired form 

error e a^Qk+i 
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This is done using the same type of calculations leading to p.23p and the following 
properties, cf. ()3.30p . 

(-29,9, + 2(1 - a)dl - ajii^ = of " 



b1 V 



Lfe+i 



\k-l 



{-2dadt + 2(1 - a)dl - dam - fl)^(log6i)'] = a)2{logb,) ^ 

We also observe that in the second sum in Lab v only Vj{l, log b, log 6i ) is important . 
The rest can be also assigned to the error. Thus matching the (1 — a) 5 like terms 
we are left with the equations 

LaVo{a,\ogb) = eo(a,log6) 

La(V;,(a, log 6, log 6i)(l - a)^) = e°(a, log fe, log6i)(l - a)^ 

Matching the (1 — a)^^ at a = 1 we get the boundary conditions (recall that 6i 
here is treated as a parameter) 

(3.36) 5b(6^-*^VS(l,log&,log6i)) = fe^-'=-ie](l, logfo, log6i), j = l,...,m 

More explicitly, p.36p means the following. Separating into monomials in log 6i we 
seek s' and {ci}f^Q so that 

4 {V-" E c, log^ b) = V-'^-^ ^ c? log^ b 

for given s and {c°}|^q. If j > A; then we set s' := s and 

(j - k)ci + {£+ l)ce+i ^4 0<i<s 



j - k 

whereas in case j = k we set s' := s + 1 and ci = -^f^ for all 1 < < s' (in 
particular, we generate extra powers of log 6 in this case and cq is not determined). 
Write 

m 

eo (a, log &) = ^ Pj (a) log^ 6 
i=o 

with Pj{a) is smooth on [0, 1], analytic close to a = 0, and Pj{a) — 0{a^). Then 
we solve the problems 

iaV^Oj =P„ j = 0,...,£(0) 

where we select a solution which is smooth at a = 1 . Using the notations of p.29p 
and variation of parameters, 

^o,i(a) = c.^o{a) + Co [ipo{a)ipi{u) - ?/'o(u)V'i(a)](pi(w))"-^Pj(w) 

J a 

where cq is an absolute, and c an arbitrary, constant. Note that around a = 0, 

^oj(a) = O(a^) + coj <po(a) 
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where Lafo ^ and Lpo{a) — a ^(1 + 0(a^)) with analytic O(a^) (as can be seen 
from a power series ansatz). Then define 

m 

Vo{a,\ogb) :=^yo,,(a) log^ 6- 

Even though this expression will in general be singular at a = 0, the singular part 
is of the form 

^(0) 

'/'o(a)^coj log-' b 

Similarly, we write 

e"j{a, log b, log 6i) = ^ ^ %,£,n(a) log^ b log" 5i 

k=Q e+n=k 

where qj.e,n are smooth, analytic around a = and vanishing to second order at 
a = 0, and solve the problems 

La[{l - a)^Vj-/,„(a)] = (1 - a)^q-,-,£,„(a) 

by variation of parameters, i.e., 

(1 - a)^Vj,e,nia) = Cj/,„ i/'i (a) 

J a 

where Cjj^n is arbitrary. Note that Vj/,„(a) is smooth around a — 1. As for the 
behavior around a ~ 0, one has 

(1 - a)^y,-/,„(a) = 0(0^) + c(po{a) 
as before. Moreover, since 

(l-a)-^i(a) = l + 0(l-a) 

we conclude that Vj_£.„(l) can be assigned arbitrary values. This is crucial with 
regard to the boundary condition (|3.36p . More precisely, setting 

(U) 

VS (a,log6,log6i) ^ F,, (a) log^ 6 log" &i 

fc=0 e+n=k 

we can satisfy the boundary condition (|3.36[) at a = 1. Generally speaking, the 
approximate solution 

m 

Vs^ng{a,b) := b-''Vo{a,\ogb) + (1 - a)^ ^ (a, log &, log6i)6J-'=6-^' 

will not be smooth at the origin a = 0, let alone vanish to fourth order. To remedy 
this problem, we subtract the correction function 

ni 

V{a,b,) b^''Vo{a,\ogbi) + (1 - a)^ ^ VS(a, log6i, log feOfer'^ 
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which solves the homogeneous equation LaV G J-k+i + and has the same 

singular behavior at a = as Vsing- More precisely, we first set 61 = & in T4ing(a, b) 
and write the resulting expression in the form 

&-'=^K(a)log''6i 
In view of our discussion regarding the singularity at a = 0, we see that 

where O(a^) is analytic and is the regular homogeneous solution, i.e., La.^\ — 0, 
(/?i(a) = a^(l + 0(a^)). Hence, we see that 

V{a, 61) hl^Y^{c^ ^o(a) + ¥^1(0)) log'' &i 
has the desired properties, i.e., 

V := Fsing - 

vanishes to fourth order at a = 0. Finally, as above one checks that 

which therefore is an error. Finally, the error /fc+i + /fc+2 generated by this entire 
procedure vanishes at least to second order at the origin as claimed. □ 

By design, Lemma 13.81 allows for arbitrary many iterations. Therefore, we can 
now carry out the process leading to vi as explained above, see p.l9|) . At each step 
we gain a power of or 6^ , while paying at most one power of log 6 and log&i. 
We iterate sufficiently often, and let 

t;2 = W2 + W3 + • ■ • 

By construction vi vanishes of order four at a = 0, therefore we can factor out an 
a* to obtain 

V2 e a4/5(l, Qi) 

Recalling also that we have neglected terms of the form (A(t)^)^^/S'*(l), we find 
that the remaining error satisfies 

as desired. 

Step 3: We now consider the general setup. Commence with e2k, k > \, 
satisfying (|3.15p and choose V2k+i so that (|3.12|) . (|3.13p hold with k replaced by 
k + 1. Note that we can move that part of e2k which belongs to 

into the next error, e2k+i- Hence we only need to deal with the part of e2k in 

IS\{\ogR)''-\Q2pk), 
which we denote as Cjj,. Proceeding as in Step 1, we then set 

{tXfv2k+iiR,a,b,bi) = e{R) / ^{R')t^elk{R' ,a,bM)R' dR' 

Jo 

-^R) J Q{R')t^elk{R' ,a,bM)R' dR' 
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Here we treat a, b, bi as constant parameters. Then it is clear that 

V2k+i e I SHR' {log R)^Q2mk+i)) 

We need to check that the error satisfies p.l3p for fc + 1 instead of k. The error 
is comprised of the terms arising from df, when one of the variables a, 6, 61 is 
differentiated, as well as the nonlinear terms. More precisely, we write 

62*:+! = N2k+l{v2k+l) + E^V2k+l + E°-V2k+1 

where the first represents nonlinear errors, the second represents dfv2k+i, and the 
third represents those constituents in 

{-d^ + d^, + -dr)v2k+l{R,a,b,bi) 

r 

in which at least one derivative falls on a, or b, or bi. It is straightforward to check 
that 

t'E'v2k+i e IS\R'{logRf, Q2/3(fe+i)) C I S*{R' {log R)'', Q'2p(k+i)) 

Next, the terms in t^E"-V2k+i are of the form 

[(1 - a^)dl + (a-i - 2a)da]v2k+i {R, a, b, 61) 

[(1 - a^)da + (a-i - 2a)]{dabidb,V2k+iiR,a,b, 61)) 

(1 - a^)tdtaRdBdaV2k+i{R,a,b, 61) - (1 - a'^)a^^RdadnV2k+i{R,a,b,bi) 

Each of these is easily seen to be in I S'^{R^ (log R)'^ , Q2/g(/c+i))- "^^"^ nonlinear errors 
are of the form 

6 2 

^(W2fe - Q^)v2k+1 + ^(3u|fc+iU2fe + 

For the term on the left, expand U2k = Q + X)i<i<2fe ^i- Using that 

v,eIS\R^Q2p), 

l<i<2k 

we check that 

t'^i^lk ~ Q^)v2k+i e I S*{R' {log R)", Q20ik+i)) 

Similarly, we get 
2+2 

^{3vl^,U2k + vik+i) e IS\R'{logR)\ Q2p(k+i)) 

Step 4 Commence with e2k-i, k > I, satisfying (j3.13p and choose V2k so that 
p.l4p . (|3.15p hold. Pick the leading order term in 62^-1, which can be written as 

k-l 

:=i?2^5,(a,&,^i)(logi?y, 

3=0 

with gj{a) G Q2/3fe- then claim that the error e\j._^ :— e2k-i — ^2k-i '^''^^ 
absorbed into e2k- Indeed, we can write 



"2^-1 



21 I /-I 2\ 1 

a e2k-i + (1 - a )e2fc_i, 
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and we have (1 — a^)Q2/3fe Q^pk- Next, rewrite 

fe-i 

t'^4k-i ^^hj{a,b,bi){logRy , hj{a,b,bi) ^ a? gj{a,b,bi) G a2Q2(fc-i)/3 
i=o 

We first seek an approximate solution 'W2k for (|3.6p of the form 

fc-i 

w^2fc ^Zj{a,b,bi){\ogRy , Zj G a'^S2(fc-i)/3 

This we then refine, iterating application of Lemma 13.81 sufficiently often to ob- 
tain V2k- To find the functions Zj we proceed inductively, starting with the largest 
power of log R. Indeed, matching corresponding powers of log i?, we get a recursive 
system. Denoting 

we calculate 

fe-i 

L°°W2k{a,b) ^Y.{ i^ogRyL'^Zj - 2{tdt)z, {tdt){\ogRy + 2{tdr)z,{tdr){\ogRy 



+ t^z, 



Z^{-d^^+dl + Ur){\0gRy] 

= Y,[{\ogRyLabZj+]{\ogRy-^Li,z,+3U - mogRy-^ LI,z,] 



k-1 



3=0 



where 



7-1 



13 13 
- 1 - - + — 



This leads to the recursive system for < j < fc — 1, 

(3.37) LabZj = hj - {j + l)Ll^Zj+i - {j + + 2)Lli,zj+2, Zk = Zk+i = 

Since hj G fi^S2(fc-i)/3 ^^'^ seek approximate solutions Zj G a^Q2(fe-i)/3i it 
suffices to take the principal part of the system (j3.37p , namely 

LabZj = hj + (j + 1)(1 + 2(a - a-^)da)zj+i - (j + l)(j + 2)(1 + a-^)zj+2 
Zk = Zk+i = 

For this we apply Lemma l3.8l to obtain approximate solutions zj G CL^Q2{k-i)i3 with 
lower order errors 

LabZj - {hj + (j + 1)(1 + 2(a - a-^)daZj+i) G a^Q2(fc-i)/3+i 

The other terms on the right hand side of (j3.37p have a similar form, 

(j + + 1 + 2(a - a-^)da)z,+i - [j + + 2)L2^z,+2 G Q'2(^k-i)p+i 
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In addition to the above error terms, by adding W2k to the approximate solution 
we have also generated errors from the nonlinear terms, which we recall are (upon 
multiplication by t^) 

where W2fc-i = Q + vi + . . . + V2k-i- We expand the first term here in the form 



-^{Q - 1 + Wi + . . . + V2k-l){Q + I+V1 + V2k-l)w2k 

with vi ~ vio + vii. First we write 



(Q + I)(Q - l)w2k = a- \ ,,^2 W2fc e IS\{\ogR)''-\ Q2kp) 



r 



2V-^ M-^ y--- I+i?2(i^)2 



£^(Q± l)i>j«;2fc, £jv^fju;2fc e /S'4((logi?)'^ \ Q2(fc-i)/3+2): 



— (3zi;2\u2fe-l + Wlk) e /5'l((l0gi?)'=-l, Q2fe/3), 



which we can absorb into e2k- The terms 

^2 

^((5± l)vioW2k 

are similar but simpler. On the other hand we recall from Step 1 that vn satisfies 
vii e IS^{R^, Q2/3+i)- Hence we obtain 

t^{Q - l)viW2k e a'lS\{\ogR)''-\ Q2(fc-i)/3+i) C IS\{\ogRf-\ Q'2ik-i)p+i). 

which we cannot absorb into e2k yet, whence we iteratively apply the preceding 
procedure to it. The remaining interactions satisfy at least 

■^{Q±l)vjW2k, 72 

and we re-iterate the preceding procedure for those which cannot yet be absorbed 
into e2k- We similarly deduce 

which can therefore be absorbed into e2k- We now re-iterate (sufficiently often) 
the procedure from the beginning of the present step for those errors which cannot 
yet be absorbed into e2fc, resulting in W2k = ^2/0' ^2fe' • ■ ■ ' '"'2fc ■ Finally, V2k '■= 
^^£q has all the desired properties. □ 

4. The analysis of the underlying strongly singular Sturm-Liouville 

operator 

In this section we develop the scattering and spectral theory of the linearized 
operator C The main tool developed in this section, which is crucial to this paper, 
is the distorted Fourier transform. The main difference between the linearized 
operator in [TT] and the one of this paper is that in [Tl] the linearized operator had 
a zero energy resonance and here zero is an eigenvalue. In both instances, though, 
there is no negative spectrum (unlike the semi-linear case [T2], where we had to 
deal with a negative eigenvalue and the resulting exponential instabilities). 

Definition 4.1. The half-line operator 

(f 15 24 



C:=- 



dR^ 4i?2 (l + i?2)2 



26 J. KRIEGER, W. SCHLAG, AND D. TATARU 

on L^(0,oo) is self-adjoint with domain 

Bom{C) - {/ e LmO,^)) : fj' e Aaoc((0, (x.)), Cf G L'((0,oo))} 

Because of the strong singularity of the potential at i? = no boundary condition 
is needed there to insure self-adjointness. Technically speaking, this means that Cq 
and C are in the limit point case at i? = 0, see Gesztesy, Zinchenko [5]. We remark 
that Co and £ are in the limit point case at i? = oo by a standard criterion (sub- 
quadratic growth of the potential). 

Lemma 4.2. The spectrum ofC is purely absolutely continuous and equals spec(£) = 
[0,(X3). 

Proof. That C has no negative spectrum follows from 

(4.1) £^0 = 0, MR)^ {iTr^Y 

with positive (by the Sturm oscillation theorem). The purely absolute continuity 
of the spectrum of £ is an immediate consequence of the fact that the potential of 
£ is integrable at infinity. □ 

We now briefly summarize the results from [5] relevant for our purposes, see 
Section 3 in their paper, in particular Example 3.10. 

Theorem 4.3. a) For each z € C there exists a fundamental system (j){R,z), 
9{R, z) for L — z which is analytic in z for each R > and has the asymptotic 
behavior 

(4.2) <l)[R,z)^Ri, 9{R,z)r^^R-i as R ^ 

In particular, their Wronskian is W{0{-, z),6(-, z)) = 1 for all z G C. We remark 
that (j>{-,z) is the Weyl-Titchmarsh solutioi^ of C — z at R = 0. By convention, 
z), 9{-, z) are real-valued for z e M. 

b) For each z G C, Imz > 0, let '0+(i?, z) denote the Weyl-Titchmarsh solution 
of C ~ z at R — oo normalized so that 

ip+{R,z) ^ z-i e'^^"' as R~^oo, Imz5>0 

If^ > 0, then the limit tp^{R, ^ + iO) exists point-wise for all R > and we denote it 
by i^+iR,^). Moreover, define tfj- {■, £_) := ?/'+(•, 0- Then i>+ {R, £,) , i^-{R,£,) form 
a fundamental system of C — S, with asymptotic behavior tp^ {R, ^ g^*?^^ as 
R ^ oo. 

c ) The spectral measure of £ is given by 

(4.3) MrfO-|l0o||2"'<5o + p(Orfe, p(e) := -Im m(e + zO)xK>o] 

TT 

with the "generalized Weyl-Titchmarsh" function 



(4.4) m(0 



■^Our </>(-, z) is the (f>(z,-) function from [5] where the analyticity is only required in a strip 
around M - but here there is no need for this restriction. 
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d) The distorted Fourier transform defined as 



for aU^>0 is a unitary operator from L'^{R+) to L'^{R+,fi) ^R®L'^{R+,p) and 
its inverse is given by 



(4.5) / ^ fiR) = fm\M2'MR) + lim / HR,OfiOpiOd^ 

Here lim refers to the L'^(R~^,fi), respectively the L^(R+), limit. 

Remark 4.4. It is best to view the distorted Fourier transform of any f G 
as a vector, namely / i-^ (g")) where a e M and g G L^(]R+,p). The inversion 
formula being 

poo 



The first term is the projection of / onto ^o, whereas the second one is the projection 
onto the orthogonal complement of 0o- We remark that 



■dR=l 



fl + i?2)4 6 



4.1. Asymptotic behavior of (f> and 9. Beginning with two explicit solutions for 
Cf = 0, namely 

, i?5 ^ , ^ -l-8i?2 + 24i?4fogi? + 8i?6 + i?8 



(l + i?2)2' ' 4i?i(l + i?2)2 



we construct power series expansions for from (j4.2p in z e C when i? > is fixed. 
A similar expansion is possible for 9{R,z). Since is it not only more complicated 
but also not needed here, we skip it. 



Proposition 4.5. For any z G C the solution <p{R^z) from Theorem \4.S\ admits an 
absolutely convergent asymptotic expansion 

oo 

^{R,z) = MR) + R'^^{R^zy^j{R^) 

The functions are holomorphic in ft = {Reu > "5} ond satisfy the bounds 

\Mu)\<^\u\'{u)-\ J>1 



J 



for all u G n. In particula^, in the region 

\cj,{R,0\^R'^MR)-Rk 

(4.6) , 

\dR^{R,0\-RH 



If a,b>0, then a <^b means that a < eb for some small constant e > 0, whereas a i<b means 
that for some constant C > one has C~^a < b < Ca 
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Proof. Write </)(i?, z) = X]^o ^''4'j{R)- The functions (pj then need to solve £(j)j = 
(pj-i- Since (po is not analytic, it is technically convenient to set 4>j{R) = R~^fj{R) 
(note that R~i is the decay of 0o)- Our system of ODEs is then, with j > I, 

C(R-lf,) = /o(i?) = 3^ 

The forward fundamental solution for C is 

H{R,R') = {MR)Oo{R') - MRyo{R))l[R>R'] 
Hence we have the iterative relation 

fj{R) = / RiiRThMRWoiR') - MRyo{R))fj-iiR')dR', 

Jo 

MR) 



(l+i?2)2 

Using the expressions for (po, 9o we rewrite this as 

fj{R)= / R^{-1-8R''^ + 24R'^logR' + 8R"^ + R'^)- 

- R'H-l SR' + 2AR^ log R + 8i?« + ^^^^ (^^^iff^^^ dR' 

We claim that all functions fj extend to even holomorphic functions in any even 
simply connected domain not containing ±i, vanishing at 0. Indeed, we now suppose 
that has these properties and we shall prove them for fj. Clearly, fj extends 
to a holomorphic function in any even simply connected domain not containing ±i 
and 0. We first show that at there is at most an isolated singularity. For this we 
consider a branch of the logarithm which is holomorphic in C \ R~ and show that 
fj{R + iO) = fj{R — iO) for R<0. Disregarding the terms not involving logarithms, 
we need to show that for any holomorphic function g we have 

j-R+iO rR-iO 

/ (logR' -log{R + iO))g{R')dR' = / (logR' -\og{R- iO))g{R')dR' 
Jo Jo 

This is obvious since for i?' < we have 

log(i?' + iO) - \og{R + iO) = log(i?' - iO) - log(i? - iQ) 

The singularity at is a removable singularity. Indeed, for R' close to wc have 
\fj-\{R')\ < \R'\ which by a crude bound on the denominator in the above integral 
leads to |/j(-R)| ^\R\ (again with R close to 0). This also shows that fj vanishes at 
(better bounds will be obtained below). The fact that fj is even is obvious if we 
substitute 2 log R' and 2 log R by log R'^ respectively log R? in the integral. This is 
allowed since due to the above discussion we can use any branch of the logarithm. 
Indeed, denoting /j_i(i?'2) = /j_i(i?') the change of variable i?'^ = v yields the 
iterative relation, with /o(u) = (i+^i 



~ r 

fj{u)= I u^{-l-8v + 12v'^\ogv + 8v^ +v^) 
(4.7) ^° 

-v^{-l-8u+ Uu^ logu + 8u^ + u^) 



2t;2(l + «)2(l + t,)2 



dv 
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Next, we obtain bounds on the functions fj. To avoid the singularity at —1 we 
restrict ourselves to the region U = {Keu > —^}. We claim that the fj satisfy the 
bound 

The kernel above can be estimated by 

- Sv + 12^2 logw + 8v^ + v^) - -8u + 12u'^logu + 8u^ + u'^) \u\ 



2t;2(l + u)2(l + t,)2 
We have 

l7,|2 

l/o(")l < 



|2 



1 + |m|2 
which yields 

l7iHI < C \u\' — ^ dx<C \u\'{u)-' 
From here on we use induction, noting that for j > 1 

\fj+i{u)\ <-^J^ x^xr'\u\'dx < ^-—^\u\^+^{u)-' 

Finally, note that the functions 0^ are given by 0j (u) = u^^ fj (u) and satisfy the 
desired pointwise bound. 

The statement (|4.6p follows from the fact that |0i(u)| > u for u ^ 1. □ 

We note that although the above series for converges for all R, z, we can only 
use it to obtain various estimates for (f> in the region 

|z|i?2 < 1. On the other hand, 
in the region ^R^ > 1 where z = ^ > 0, we will represent (p in terms of tp^ and use 
the tp^ asymptotic expansion, described in what follows. 

4.2. The asymptotic behavior of ip^. The following result provides good asymp- 
totics for ip^ in the region R^(_ > 1. 

Proposition 4.6. For any £, > 0, the solution ip^{-,£,) from Theorem \4-S\ is of the 
form 

^+{R,Q=^-ie'^^^aiR^^,R), i?^^ > 1 
where a admits the asymptotic series approximation 

°° 15i 1 

3=0 

with zero order symbols (R) that are analytic at infinity, 



^3 

film IC/?/9o^'=7^^ 



sup|(i?5K)V+(^)l <oo 



R>0 

in the sense that for all large integers jo, and all indices a, pi, we have 

jo 

< r a - n^^°~^ 



sup 

R>0 



3=0 



for all q > 1. 
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Proof. With the notation 

we need to solve the conjugated equation 

(4-8) {-dl - + ^ - J^^) -m^R) = 

We look for a formal power series solving this equation, i.e., 

oo 

(4.9) a(q,i?)=^r*/,(i?) 
This yields a recurrence relation for the /j 's, 

m^) jTTiw) /.-(^)' /" = 1 

which is solved by 

Extending this into the complex domain, it is easy to see that the functions fj are 
holomorphic in C \ [—«,«]. They are also holomorphic at cxd, and the leading term 
in the Taylor series at oo is R^^ . At one has the estimate 

\{Rdn:ffj{R)\<c,kR~' yR>0 
which is easy to establish inductively. The functions 

now satisfy the desired bounds due to the bounds above on fj. The remainder of 
the proof is the same as in our wave-map paper [11] and we skip it. □ 

4.3. The spectral measure. We now describe the spectral measure by means 
of (|4.4p . This requires relating the fmictions (f>, 9 and ip^. By examining the 
asymptotics at i? = we see that 

(4.10) W{e,(j))^l 

Also by examining the asymptotics as i? — > oo we obtain 

(4.11) W{il)+,i}~)^-2i 
Lemma 4.7. a) We have 

(4.12) </)(i?, = a{S,)^+{R, + a(C)^+(i?,e) 
where a is smooth, always nonzero, and has size 

KOI 1^-1 ,^^>i 

Moreover, it satisfies the symbol type bounds 

\m)'a{0\<ck\am ve>o 

b) The absolutely continuous part of the spectral measure /^(rf^) has density p(^) 
which satisfies 

... / 1 e « 1 
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with symbol type estimates on the derivatives. 

Proof, a) Since (f> is real- valued, due to (|4.11[) . the relation (|4.12p above holds with 

We evaluate the Wronskian in the region where both the ip^{R,^) and ^) 
asymptotics are useful, i.e., where R^^ « 1. The bounds from above on a and its 
derivatives thus follow from Propositions 14.51 and [ 
For the bound from below on a we use that 



2\dR^+{R,0\ 



which was obtained in [TT] . We use this relation for R ^ 2 with a small constant 
i5. Then by Proposition 14.51 we have 



\dR<l>{R,0\ > 

while by Proposition [46l 




This give the desired bound from below on a. 
b) In [11] it was shown that 

TT 

The bounds on p(^) now follow from part a). □ 



5. The transference identity 

We now write the radiation part e in terms of the generalized Fourier basis 
(j){R,£^) from Theorem 14.31 i.e.. 



eiT,R)^xo{T)MR)+ / x{T,0<P{R,OpiOd^ 
Jo 

As in [11] , [12] we define the error operator /C by 
(5.1) Rd^ = -2id(u + Ku 

where the hat denotes the "distorted Fourier transform" and the operator —2^9^ 
acts only on the continuous part of the spectrum. In view of Remark l4.4l we obtain 
a matrix representation for IC, namely 



K = 



c. 
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Here 'c' and 'e' stand for "continuous" and "eigenvalue" , respectively. Using the 
expressions for the direct and inverse Fourier transform in Theorem 14.31 we obtain 







MR), MR)) , 


ICecf = 




f{ORdRHR,OpiOd^,MR)),, 


ICceiv) = 


(ror^ 


MR)A{R,v)) , 

' ^R 




(i: 


f{ORdRm,Op{i)di,^{R,v)) , 

' ^R 



(5.2) +(/ 2^d^f[0<t>{R.0p{i)d^,m,v)) , 

Wo ' i-R 

Integrating by parts with respect to R in the first two relations we obtain 

1 1 r°° 



where 

KM^{RdRMR),<t>{R,Ti), , 

Integrating by parts with respect to ^ in (j5.2p yields 



(5.3) 



JCccf{v) = { / fmRdR-2^d^]<l>[R,^)p[0d^,<l>[R,v)]^^ 



-2(1. ^)/(,) 



where the scalar product is to be interpreted in the principal value sense with 

/eCo-((0,oo)). 

In this section, we study the boundedness properties of the operator JC. We 
begin with a description of the function Ke and of the kernel Ko{r], ^) of JCcc- 

Theorem 5.1. a) The operator JCcc can be written as 

(5.4) /C,, = -(^ + ^)<5(^-r;) + /Co 
where the operator ICq has a kernel Ko{ri,^) of the /orrrQ 

(5.5) Ko{v,0 = -^m,v) 

with a symmetric function F{^, rf) of class in (0, oo) x (0, oo) and continuous 
on [0, oo)^. Moreover, it satisfying the bounds 

1 e + r; < 1 

(e + 77)^^(1 + le^ e + ??>i 



< 



^The kernel below is interpreted in the principal value sense 
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where N is an arbitrary large integer. 

b) The function Ke and K'^ are bounded, continuous, and rapidly decaying at 
infinity. 

Proof. We first establish the ofF-diagonal behavior of /Ccc, and later return to the 
issue of identifying the (5-nicasure that sits on the diagonal. We begin with (|5.3p 
with / e C^((0,oo)). The integral 

Jo 

behaves like at and is a Schwartz function at infinity. The second factor 
<j){R,ri) in (|5.3p also decays like i?2 at but at infinity it is only bounded with 
bounded derivatives. Then the following integration by parts is justified: 



Moreover, 



poo pOC 

{Cu){R)= fmC,RdR]cj){R,Op{Od^+ f{0{RdR-2^d^mR,0pi0d^ 

JQ Jo 

OO pOG 

f{0[C,RdR]cl,{R,Op{Od^+ / ^f{0(RdR-2^d^)c^{R,0p{0d^ 
Jo 



-2 ^fimR.OpiOd^ 

Jo 

with the commutator 



48 Q67?^ 

[A Rdn] = 2C + - ^^^^ =: 2C + UiR) 

Thus, 

pOO pOO 

{Cu){R)^ fm{R)HR,Op{Odi+ ^fmRdR-2^d^)ct>{R,0p{0di 
Jo Jo 

Hence we obtain 



VlCccfiv) - iCccmiv) = { / fmiR)^{R, OpiO d^, m v) , 

The double integral on the right-hand side is absolutely convergent, therefore we 
can change the order of integration to obtain 



iv ~ OKoiv,^ = p{0{u{R)<PiR,0,<f>iR,v) . 

This leads to the representation in (j5.5p when ^ 7^ // with 

Fitv) = {u{RmR,0,(biR,v)] 

It remains to study its size and regularity. First, due to our pointwise bound from 
the previous section, 

<min(i?i(i?)-4(l + i?4^),e-3) VO<C< 1 

10(^,01 <min(i?^r^) V^>1 
Note that these bounds imply that for all ^ > 0, 

(i?)-2|0(i?,C)|<0o(i?)<(i?)-^ 
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Hence, \F{£^, ri)\ <1 for all < ^, 77 < 1. Moreover, F{^, rj) is continuous on [0, 00)^ 
by dominated convergence. Finally, using that \<j){R,^)\ < when ^ > 1 implies 
that 

(5.6) \F{^,v)\ < {0~Hv)--' VC,??>0 

We shall improve on this in a number of ways, but first we consider derivatives. By 
the previous section. 

Consequently, if < ^,7/ < 1, then 

/•OO 

Jo 

_ 1 

PT] poo 

< {R)-^{\ + R'^ri)dR+ ^{R)-^-i dR<l 

Jo Jrj^i 

whereas ifO<^<l<r7, then 

/•OO 

\deF{^,r])\< {Ryirj-i dR<-n~i 
Jo 

If < ?7 < 1 < ^, then 

\d^F{^:V)\ < / {R)~'^mm{Ri,R^-i)mm{{R)-^l + R^r]),7]-i)dR 
Jo 

_ 1 

<r^ / {RyHi + R'^r])dR+ _yRy^Rr'^v^^dR<c^ 

Jo Jr]^2 

Finally, for 1 < ^,77, 

/•OO 

Jo 

To summarize, 

(5.7) \d^Fi^,v)\<{0-Hv)-K \dnF{^,v)\<{0-Hv)-'^ V^,77>0 
For the second derivatives we use that 

\dl(l,iR,0\<mm{R'^,R^^-i) Ve>l 

\dl(l){R,(,)\<mm{Ri,R^r^) V < < 1 
which imply the bounds we always have the estimates 

\dlr,F{^,v)\<riv~-^ Ve>l, 77>1 

(5.8) |5|F(^,77)| <r^r;-t Ve>l, 77>1 

The bounds (|5.6p . (|5.7p . and (|5.8p are only useful when ^ and 77 are very close. To 
improve on them, we consider two cases: 
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Case 1: 1 < C + '7- To capture the cancellations when ^ and 77 are separated we 
resort to another integration by parts, 
(5.9) 

rjFitr]) = (uiR)<j>iR,0,^HR,v)) = ([C,U{RmR,0,(biR,v)) + ^Fi^,Tj) 



Hence, evaluating the commutator, 

(5.10) (77 - OFit V) = -(i2UR{R)dR + URR{R))(f>iR, 0, ^R, v)) 

Since Ur{0) = it follows that {2UR{R)dR + URR{Rj)(j){R, ^) vanishes at the same 
rate as '/>(i?, f ) at i? = 0. Then we can repeat the argument above to obtain 

(77 - O^FiC, v) ^~{[C., 2URdR + Urr]^{R, 0, 0(i?, v)) 
The second commutator has the form, with V{R) := —24(1 + R^)^^, 

[£, 2URdR + Urr] = AUrrC - AUrrrOr ~ Urrrr ~ 2U rVr - AUrrV 
Since V{R), U{R) are even, this leads to 

iv - O^Fi^, r,) = ((JJ'^''\R)dR + C/™^"(i?) + eC/™^"(i?))</'(i?, 0.4>{R. V)) 

where by 1/°'^''', respectively [/'^■"^"^ we have generically denoted odd, respectively 
even, nonsingular rational functions with good decay at infinity. Inductively, one 
now verifies the identity 

(5.11) 

k-l k 

iv - e'Fi^, '7) = ( ( E + E ^'urr (r)) m o, m v)] 

j=0 1=0 

{R)Mf{R)\ + \u^,riR)\ < {R)-'-^' V J, ^ 

By means of the pointwise bounds on (f> from above as well as 

max((i?)-5,^3) < 1 if 0<C<1 
mm{Ri,^-i) < ^"i if C > 1 



we infer from this that 



\FiLv)\<^\^§r- VC, r7>0 



Combining this estimate with (j5.6|) yields, for arbitrary N, 

(5.12) \F{tv)\<i^ + vrH^ + \^^' provided £ + r7>l, 

as claimed. For the derivatives of F we follow a similar procedure. If £ and 77 are 
comparable, then from (|5.7p . \djjF{^, ri)\ < (C) We will use this bound only when 
1^5 — 775] < 1 which of course implies that £, ^ rj >1. Thus, we now assume that 
1^2 — 7y 2 1 > 1 which is the same as |C — J?! > {(. +v)^ ■ 1^. this case, we differentiate 
with respect to rj in (I5.11|) . This yields 

k-l k 

(77 - O^'d.Fi^, 77) = (( E e Kf'iR) dn + Y. ^'UtriR)) <P{R, 0,d,HR, r;)) 
j=o e=o 

^2k{7j-0"'-'F{tri) 
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Using the bound on F from Ij5.12p as well as the usual estimate on drj(l){R, rj), leads 
to 



(5.13) 



-N 



provided S, + Tj > 1 



The second order derivatives with respect to ^ and rj are treated in an analogous 
manner. We note that it is important here that the decay of Ul^j'^ and J/^l"^" 
improves with k. This is because the optimal second derivative bound for small 77, 
viz. \djj(j){R, rj)] < R2 , has a sizeable growth in R. 

Case 2: ^,77 <C 1. First, we note that 

F{0, 0) = {U^o, 0o) = {{[C, Rdi^] - 2/:)<^o, 0o) = 
Together with the derivative bound (|5.7p , this implies that 

as claimed. To bound the second order derivatives of F we recall the pointwise 
bounds, for < ^ < 1, 

\di(f>{R,0\ <mm{R^,R^-i) 
IfO<i^<77<l, then these bounds imply that 
(5.14) 



\doim,v)\ 



< 



e"5 



{R)-^R''dR+ / ^ {R)-^R^r]--i dR+ / ^{Ry\ir])-^ dR 



i 5 



This bound is only acceptable as long as ^ and 77 are comparable. Otherwise, if 
< ^ ^ 77 < 1, then one needs to exploit the oscillations of (3^0(i?, 77) in the regime 
-R^77 > 1 as provided by Proposition 14.61 and Lemma 14.71 Thus, write 

dr^m.ri) = 9„[a(77)7/'+(i?,7/)+a(77)7A+(i?,7,)] = 2 Re a4a(7;)77-^ e^^"' (T(i?77^ i?)] 

= 2 Re [(0(77)77-3 )'e'-R'?* a(i?77^ , i?)] + i? Re [«a(77)77-ie*-^''* cr(i?7/5 , R)] 

+ i?Re [a(77)77-3e''-«"*cr,(i?775,i?)] 



Therefore, 

(5.15) 
(5.16) 
(5.17) 



^ UiR)d^q^iR,Od^q^iR,v)dR 



< 



^ U{R)d^<t>{R,i){a{r,)r, 



cr(i?775,i?) dR 



^ 1 

^ RU{R)d^(t)iR,Oair])T]-i a{Rr]i , R)dRe'^'^^ dR 

° - 1 J- 

^ RU{R)d^(f){R,0aiv)v^^crg{Rr]2,R)dRe'^'^^ dR 



The term on the right-hand side of (15.15^ is bounded by 



^ R ^rj dR<ri 
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whereas (|5.16|) and ()5.17|) require integrating by parts. It will suffice to consider 
the former. Using that \dii(,(f>{R,£,)\ ^ min (i?3,i?^-4) and that \dqa{q,R)\ < R-\ 
we obtain 



J, 2 



^ 1 

^ RU{R)d^(l){R,C)a{'q)r]-ia{RT^'^,R)dR,e'^'^^ dR 



< 



V-^\RU{R)d^<l>{R,0\ 



+ , [{R)-'\d^HR,0\ + {Rr'\dRiHR,0\]dR<v-' 

In conclusion, for all < ^, ?/ < 1, 

as desired. Next, consider d^F(^,r]). The bound 

\dlF{tv)\< / {R)~^mm{R^,R^r^){R)^ dR<C^ 
Ja 

is acceptable as long as < < ^ < 1. If, on the other hand, < ^ ^ 77 < 1, then 
differentiating in (|5.10p we obtain 



(77 - 77) = 2d^Fit V) - 0, {2URdR + URR)(b{R, 77); 

which implies that 

7j)\ < 77-1 [\d^F{t 77)1 + 1 0, URRcj^iR, Tj))\ + \{dlcj,{R, 0,R-'URRdRC^{R, 77) 

The first term in brackets is ^ 1, the second is bounded by 

_ 1 

/ {R)-^R)^R)-i{l + R'^r])dR+ ^R-^R^~idR<l 

Jo Jri^^ 

whereas the third is the same as the second in the range < i? < 77^^, whereas in 
the range R> rj^2 we need to integrate by parts; schematically, this amounts to 



< 1 



The full details are essentially the same as in the previous integration by parts step 
and we skip them. 

Next, we extract the 5 measure that sits on the diagonal of the kernel K from 
the representation formula (|5.3|) . see also (|5.4p . To do so, we can restrict ^,77 to 
a compact subset of (0, 00). This is convenient, as we then have the following 
asymptotics of for i?^2 ^ l: 



(t){R,Cj = Re 



0{R- 



{ROr - 2^d^)(l>{R, = -2 Re 



+ 0(i?-2) 



where the O(-) terms depend on the choice of the compact subset. The R ^ terms 
are integrable so they contribute a bounded kernel to the inner product in ([57 
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The same applies to the contribution of a bounded R region. Using the above ex- 
pansions, we conclude that the (5-nieasure contribution of the inner product in ()5.3|) 
can only come from one of the following integrals: 

(5.18) 



fiOxiR) Re 



/o ^0 
(5.19) 
1 

^ 2 
(5.20) 
1 



^a,(a(Or ^)a(^)^-^e«(«*+"*) (l + (l + 



oo /'OO 



fiOx(.R) e55(a(Or^)«(r/)ry-^e'^(«'-''') (l + 



15i 



1 - 



8Rf] 
15i 



p(0 d^dR 



8Rr]^^ 



p(0 d^dR 



I r r/(0x(i?)C9e(a(e)r^)a(^)^-^e-«(«*-''^)(l--^)(l + -^) piOd^dR 
^ Jo Jo ^ 8it^2/\ 8it772 / 

where x is a- smooth cutoff function which equals near R = and 1 near R = oo. 
In all of the above integrals we can argue as in the proof of the classical Fourier 
inversion formula to change the order of integration. Integrating by parts in the 
first integral (|5.18p reveals that it cannot contribute a ^-measure. Discarding the 
R~^ terms from (|5.19p and (|5.20p reduces us further to the expressions 



(5.21) 



OO /'OO 



/(e)x(i?) Re 



p(^) d^dR 



(5.22) 
15 



OO poo 



/(e)x(i?)Ini 



Caaa(e)r^)a(r?)r?-^e'^(«"-''*) R-'iC'^-V-hpiOd^dR 



The second integral (|5.22p has both an R^^ and a (C~^ — '7"^) factor so its contribu- 
tion to K is bounded. The first integral (|5.2ip contributes both a Hilbert transform 
type kernel as well as a (5-measure to K. By inspection, the S contribution is 



1 
2 

— TT Re 



Re 



piOdR 



^2^C3p(e)Re e9c(a(e)r^)fl(0r' S{^~7j) 



2<3p(^)Re 



--r^|a(0P+^^«(0«'(e) 



p(0 



<5(e - r?) 



where we used that p{0^^ = M'^l'^ the final step. Combining this with the 
(5-measure in (|5.3|) yields (|5.4[) . 

b) Arguing as in part (a) we have 



Keiv) 
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For F we use the representation in (I5.11|) with ^ replaced by and (f>{-,£,) replaced 
by 4>Q. The conclusion easily follows from pointwise bounds on (j){-,ri) and its 
derivatives. □ 

Next we consider the mapping properties for /C. We introduce the weighted 
spaces L^'" of fimctions on spec(£) with norm 

/•OO 

(5.23) ll/ll^,»:=l/(0)P+ / l/(C)P(0'>(Orfe 

" Jo 

Then we have 

Proposition 5.2. a) The operators ICq, K. map 

ICq : Lp' ^ Lp' ^ ' , K. : L^' . 

b) In addition, we have the commutator bound 

with acting only on the continuous spectrum. Both statements hold for all 

Proof. We commence with the /Co part, a) The first property is equivalent to 
showing that the kernel 

pHv){vr^'^'Ko{v..O{0-''p-^iO ■■ i'(M+) ^ L2(R+) 
With the notation of the previous theorem, the kernel on the left-hand side is 

We first separate the diagonal and off-diagonal behavior of Kq, considering several 
cases. 

Case 1: (C,??) G Q [0,4] x [0,4]. 

We cover the unit interval with dyadic subintervals Ij — [2^^^, 2^+^]. We cover 
the diagonal with the union of squares 

2 

A= [j I,xl, 

j=-oo 

and divide the kernel Kq into 

^qKq = ^AnqKo + ^q\aKo 

Case 1(a): Here we show that the diagonal part lAnqKo of Kq maps to L^. 
By orthogonality it suffices to restrict ourselves to a single square Ij x Ij . We recall 
the Tl theorem for Calderon-Zygmund operators, see page 293 in [16 : suppose the 
kernel K{r], ^) on defines an operator T : S ^ S' and has the following pointwise 
properties with some 7 G (0, 1] and a constant Cq: 

(i) \K{r,,0\<CQ\^-v\-' 

(ii) \K{r^,0 ~ K{r^',0\ < Cq\v - VI^IC - v\-'-^ for aU [77 - < \^ - rj\/2 

(iii) \K{r^,0-K{r^,C)\<CQ\^-C\^\^-v\-'-'' i0Te.\l\C-C\ < 1^-^1/2 
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If in addition T has the restricted boundedness property, i.e., for all r > 
and ^0,% e R, ||r(w'^'«")ll2 < C^r^ and ||r*(w'-^"«)ll2 < Cqt^ where = 
— Co)/'') with a fixed bump-function w, then T and T* are i^(R) bounded 
with an operator norm that only depends on Cq. 

Within the square Ij x Ij , Theorem 15.11 shows that the kernel of Kq satisfies 
these properties with 7 = 1, and is thus bounded on . 

Case 1(b): Consider now the off-diagonal part Igy^ifo. In this region, by 
Theorem 15.11 

\K^{ri,i)\<{iri)--- 
which is a Hilbert-Schmidt kernel on Q and thus bounded. 

Case 2: (C,?/) € Q^. We cover the diagonal with the imion of squares 



and divide the kernel K^) into 



Case 2a: Here we consider the estimate on _B. As in case la) above, we use 
Calderon-Zygmund theory. Evidently, \Ko{r],£^)\ < 1^ — r]\^^ on B by Theorem l5.1l 
To check (ii) and (iii), we differentiate Kq. It will suffice to consider the case where 
the i9{ derivative falls on F{(^,r]). We distinguish two cases: if 1^3 — | < 1^ then 
[■f ~ ^1 ^ which implies that 
f~5|£-£'| 

I J I ^ TT-^ y\^-e\<\^-v\/2 

if, on the other hand, j^^ — 7y^| > 1, then 

,/'''5r^'' < V le - ri < - 

which proves property (iii) on B with 7 = ^, and by symmetry also (ii). The 
restricted property follows form the cancellation in the kernel and the previous 
bounds on the kernel. Hence, Kq is bounded on B. 

Case 2b: Finally, in the exterior region \ i? we have the bound, with arbi- 
trarily large N, 

|i^o(^,e)l< (1 + 0^^(1 + ^)-^ 
which is bounded by Schur's lemma. 

This concludes the proof of the first mapping property in part (a). The second 
one follows in a straightforward manner since K,, is rapidly decaying at oo. 

b) A direct computation shows that the kernel Kq°"^ of the commutator [^d^ , Ko] 
is given by 

interpreted in the principal value sense and with i^^"™ given by 
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By Theorem 15.11 this satisfies the same pointwise off-diagonal bounds as F. Near 
the diagonal the bounds for and its derivatives are worse than those for F by 

a factor of (1 -l-^)^ . Then the proof of the commutator bound for Kq is similar 
to the argument in part (a). 

The remaining part of the commutator [/C,^9^] involves 

(i) The commutator of the diagonal part of ICcc with ^9^. This is the multipli- 
cation operator by 

which is bounded since p has symbol like behavior both at and at oo. 

(ii) The operator ^d^K-ce which is given by the bounded rapidly decreasing func- 
tion ^d^K,{^). 

(iii) The operator ICec^d^ given by 

Jo Jo 
which is also bounded due to the properties of Ke- □ 

6. The second order transport equation 

This section is devoted to the study of the linear problem p.2p which we restate 
here in the form 

(6.1) {~df + df. + r-'dr + 2r-\l 3Q{Rf))e = f- 12u;'^^Si^_^e 

We recall that the second term on the right-hand side here arises due to fact that 
its decay is of the same nature (namely w^) as that of other error terms which we 
will encounter in the parametrix construction of this section. By doing this, the 
remaining terms in the nonlinearity Af in (j2.2[) decay more rapidly at infinity. Our 
main result asserts that 

Proposition 6.1. The backward solution e for (|6.ip satisfies the bound 

(6.2) \\e\\Hl,<^\\f\\Ll 
for all large enough N . 

Proof. We work in the coordinates {R, r) given by 

R^rX{t), T = ^ X{s) ds = if] + l)-^\\ogtf+^ 

for any < < < 1. For future reference, we note that 

tX{t) = ((/3 + l)r)'^, A(r) = ((/? + l)r)7'^e»''+i)^)^ 
We introduce the auxiliary weight function uj(t) 

u;ir) := X-'X^ir) = j^r'^ + HP + l)r)-^ 

and note that 

P -1 



(6.3) (iA)-i = lu{t) 



f3 + l 
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Then 

dt = ^{dr + Rrdn) - -A(r)(5, + ujRdR) 

= A2(t) [{dr + LoRdnf + Lj{dr + LjEOr)] 

therefore the equation (|6.ip becomes 

^ /-l^'^ (l + i?2)3 ^ 

At this point it is convenient to switch to the notations 

(6.4) £(r, R) = R^e{T, R), J{t, R) R^\-^![t, R) 

Since 

+ ujRdR)R-^ =dr+ ujRBr - cj/2, 
one concludes from conjugating the previous PDE by i?2 that 



15 



24 



(6.5) 

where cu := dr^ and 

^ ' 4i?2 (1 + i?2)2 

Written in terms of (e, /) the estimate (j6.2p takes the fornfl 
(6.6) 
where 



R^(l~R^). 



(l + i?2)3 



lell^, = sup r^-i-7^||£(r)|U. +T^-i(||/:^e(r)|U. + ||(a.+c.i?ai^)£(r)|U.) 



-'•>To 



respectively 



sup r^||/(r)||i. 

r>To 



In order to take advantage of the spectral properties of the operator C we conju- 
gate the equation (j6.5p by the Fourier transform J- adapted to C The transference 
identity is 

TROrT-^ = -2^8^ +K, 

where 



cc 



:Id 



/Cec 

/C,e -(l + 7yp'(?7)Mr7))(5(C-r])+/Co 



and 



/•OO 

Jo 



^Here we slightly abuse notations since the A^'s in 116.21 1 and 116.61 1 do not coincide, instead they 
are linearly related. 
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We write 



where 





2e9c + (l + ep'(e)M0). 

ICec 
ICce ICo 








0" 







So 



Then 



T{dr + UjRdR)T-^ =Dr + -- iolCnd, = dr - Lo{l + JCd) 



therefore 

T{dr+UjRdRfT-^ = {Dr + 2ujlCndDr+ Uj''{[JCndXd] + - lCnd)~ ^ICnd 

Next we consider the Fourier transform of the last term in (|6.5p . which we express 
in the form 



ce \J cc 



We note that 



_2,, -12i?2(l_i?2 



(l+i?2)3 



1\ 1 3 



6 / 10 5 



while 



and 



— 12i?'^(l — i?2) 



Jo 

We remark that the kernel Je is bounded and rapidly decreasing at infinity. Finally, 

POO 



with 



12S^i-^j5p(r7)0(i?, ?7)0(i?, OdR 



(l + i?2)3 

This is bounded and has the ofF-diagonal decay property 

(6.7) Jcc(e,r?)|<(l + 0"^(l + ie^-^^l)-'^ 

Taking into account all the notations above, the equation (|6.5|1 becomes 



Te = TJ — 2ujJCndDrJ^£ 



Next, write Te 



and Tf = 



, or equivalently. 



i{T,R) = xo{t)MR)\\M\2^ + / x{T,0^{R,Op{OdC=-£o+Sc 
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where Iq _L Ic for all r > (recall (?!)o(i?) = i?5 (1 + i?2)-2 ^ ^^-^.g 

a; 

2 



the system for 



we compute 



K^dK^nd 

We also note that 



JC 

K-ce ICq 

ICec 
K-ce K-o 



Bo 



1 ^ ro -JCecBo 

SoJ [O -/Co6o 

JCec] _ r 

JCce Kq —BqICcc —BqICq 



-K K 



WMh 



6: 



17 13 

120 ^ 4 ~ 5 



Then we seek to write the equations for xq and x in the form of a diagonal system 
with perturbative coupling, 



p 


Xq 




90 


+ Q 


xo 




X 




9 




X 



Q = 



■ 

Qce Qc 



(6.8) 
where 

with the principal part given by 

Pe = -dr{dr - Uj) 

respectively 

Pc = -Dl -uDr-i 

and the coupling terms of the form 

QecX = Up'TZecX — 2ujlCecDrX 

with 

Tlec = [lO~'^UJ — 1)/Cec + Kec^o + Jec 

while 

QceXo = OJ^TlceXo — 2ujlCcedTXo, QccX = Up'TlccX — 2ujlCoDrX 

with 

Ti-cc — [ICo, Bo] + /Cq + ICcelCec + Jcc 

respectively 

Our main solvability result in Proposition 16.11 for the equation (|6.ip is restated 
in terms of the system ()6.8p as follows: 

Proposition 6.2. For each with (90,9) which satisfy 

\9o{r)\<r-^, \\9ir,-)\\Ll<r-^ 

there exists an unique solution (x, xq) for the system (|6.8p decaying at infinity. This 
solution satisfies the bounds 



1 



2/3+1 
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respectively 

(6.10) ||x(r)|U. < lr-^+^, ||e^a;(r)|U. + p.x(r)|U. < 1 

Proof. Our strategy is to solve first the simpler linear equations 

(6.11) -dr{dr - uj)xo = go 



.12) 



Dl - ujDr - e 



Then we will show that the right hand side in the system (|6.8p is perturbative. We 
start with the linear operator governing xq , and introduce the appropriate function 
spaces for xq and go: 

II5o||y« = sup r^|.go(T)| 

T>To 

WxoWx- = sup r^-^|a;o(T)|+r^-7^|9.Xo(T)| 

T>Ta 

Lemma 6.3. The backward solution operator xo — T^go for (jO.lip satisfies the 
estimate 

(6.13) WTegoWx- < \\9oi<j)\\y. 

for any N > 2. 

Proof. A fundamental basis of solutions of —dridr ~ lo) is given by 

/>OC 



and has Wronskian W{t) — A(t). Then the backward fundamental solution is given 

by 

(6.14) Uo{T,a)^W-^{a){a+{T)a^{a)-a+{a)a^{T))^\{T) j \~\s)ds 
A direct computation shows that Uq satisfies the bounds 

\Uoir,a)\<4^, |a.C/o(r,a)|<r-^ + '^(")"(") 



uj{t) ' ' ^ 'J-'(o') A((t) 

The conclusion of the lemma follows. 

□ 

Next we bound the solution of the equation (|6.12p . which is hyperbolic. One is 
tempted to define spaces and in a manner which is similar to and . 
This would work for the linear theory for (|6.12p . but would not be strong enough 
in order to treat the right hand side in (|6.8p in a perturbative manner. Instead we 
define some stronger spaces using the additional weight 

m(o-r+r'^ 

where > is a fixed small parameter. We define the space L^L^ with norm 



and the dyadic L space l^Lp^ with norm 



= sup ||cr^ "<'' + ^'g(cr)||L2 ([^,2r]xl 



T>To 
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Then we define the Y space as a sum of two spaces, 
\\g\\Y« = inf \\gi\\L^Ll + \\g2icr)\\ 

9=91+92 " P 

Similarly we introduce the space with norm 



,2 + \\{^2x,D^x)\\l^^ 



Then our solvability result for (|6.12p is as follows: 

Lemma 6.4. The backward solution operator x ~ Teg for the equation ()6.12|) 
satisfies 

(6.15) ||Te.g||x« < My- 

In addition we have the smallness relation 

1 



(6.16) 
for large N . 

Proof. The equation (|6.12p is equivalent to 



(6.17) 



We substitute the functions {x,g) by {y,h) where y — p^x and h = p'^g. This has 
the effect of removing the weight p from the estimates. The functions {y, h) solve 

2 



(6.18) 



(dr ~ 2uj^d^ 



The characteristics of the homogeneous operator on the left are (r, A ^(t)^o) which 
means that 

{dr - 2u;^d^) f{r,0 - ^fir,ar)), C(r) = X-'{r)^o 
Hence, we are reduced to solving the ODE 



(6.19) 



d-^ + LoiT)dr~X-'{T)^o A-2y(T,e(T)) =A-2/j(r,e(r)) 



with ^0 > fixed. The homogeneous equation has exact solutions 



This is no surprise since the equation (|6.12p is equivalent to the constant coefficient 
wave equation in the t, r coordinates. 
Since the Wronskian 

it follows that the backward solution to (|6.18p has the form 

y{r,Co)^^o'J ^sin(eo'y X-\u) du) hia,^^)) 
Define the forward Green function 

U{t, a; sin (^U(r) £ X-\u) du) 



da 
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Since go = ^'^i''')^, ^(0") — C^^i^)^ ^(o"); we can write 



U{T,a;Oh{<j,a<j))da 



To estimate D-rH it is also convenient to evaluate directly 

(6.20) DrU {t, a; = cos (^^^X{t) J \-^{u)du^ 

To estimate the solution y we either bound | sin(w)| < \v\ or | sin(ti)| < 1. Using 
that 



A(r) / X-\u)du<uj-^{T) 



one obtains 
(6.21) 

as well as 
(6.22) 

We denote 
(6.23) 



A(cr) 



C'\UiT,a;0\ + \DrUiT,a;0\ 



< A(r) 



An immediate consequence of f|6.2ip and (|6.22p is the estimate 

X-\<j)\h{a,aa))\da 
From this we need to conclude that the following four bounds hold: 



(6.25) 

(6.26) 
(6.27) 

respectively 
(6.28) 



1 



Z 100 7-2 



Zll/oc 7-2 
'"jv-i-'^i/ 



N' 
< J-l 

< \\h\\i- 



< 



\h\l 



Taking norms in ^ on both sides of (|6.24p we obtain 



Mr)\\ 



< 



\h{a)\\i^2 da 



which leads directly to (|6.25p . 

Adding flow invariant weights to the above bounds we get 

-1 f^Kr) 



hia] 



and by Cauchy-Schwarz 



N 



00 



A(cr 

-1 f^Kr) 
X{a) 



da 



L2 



hia] 



da 
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Hence 

f-2Ti 

r^(^-i)-^||z(r)||i. dr 



<-y a^min{a,2ri}^-i|lMa)|li.da 



where 

roc ^ pOO 

M = sup/ T~/J+TTO"^(^A(r))dT w sup / m~^(^A)— w 1 



This concludes the proof of (|6.26p . 

We now turn our attention to (|6.27p . for which we need to take h G i L. 



2 



T+75 



From (|6.24p by Cauchy-Schwarz we obtain 

X-'ia)m{a'j))oj'\<T)\hia,a<^))\da- J u:{a)m~\£,{cT)) da 
The second integral has size 0(1), therefore 

\-\<j)m{£,{a))u-\a)\h{<j,^{a))\d<j 
Hence integrating with respect to f to obtain 



OO nOO 



Mr)\\h<J miOu;-\a)\h{a,0\'d^da 
This directly implies that 

which gives (|6.27|) . 

Finally, for (I6.28p . from (|6.29p by using again Cauchy-Schwarz we obtain 

> poo 

^(r)A-2(r)|z(r,e(r))pdr< / A-2(a)m(e(a) da 

Jo 

and integrating with respect to 

) /"OC 

^{r)\\z{r)\\l. dr< / cj-\a)\\h{a)\\l. da 

Since the equation (|6.19p is solved backward in r, we can add any nondecreasing 
weight in the above estimate. In particular we obtain 

/ ^2(Ar-i)-^||^(^)||2^ / min{a,2Ti}2(A'-i)+FfT||/,(^)||2^d^ 

J Ti ^^^^ J Ti ^ 

Hence (jOS]) follows. □ 
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It remains to show that the right hand side terms in (16. 8|) are perturbative. We 
solve the equation (|6.8p iteratively and seek a solution as the sum of the series 



(6.30) 



(TQ)'^ T 



It remains to establish the convergence of the above series. By Lemmas 16.316.41 the 
backward solution operator T for P, given by 

'Te 0' 
T, 



T = 



is bounded 



T : F(f X y 



N 



X X 



N 



Hence an easy way to establish the convergence of the series in (|6.30p would be to 
show that 

WQWx^ xX^ ~.Y^ xY^ < 1 

We can establish such a bound for certain components of Q, but as a whole Q is 
not even bounded in the above setting. Lacking this, a weaker but still sufficient 
alternative would be to prove that 



\\TQ\\ 



< 1 



This is still not true, but we will establish a weaker bound, namely 
(6-31) ll^<9llxo"xx"-->Xo"xAr" ^ 1 

This ensures that all the terms in the series in (jOn)) belong to X^ x X^ . In order 
to ensure convergence we will split Q into two parts, 

Q = Qg+Qb 

The good component Qg contains most of Q and satisfies a favorable bound 



(6.32) 



< 



1 



5>Q 



Here the constant on the right can be made arbitrarily small by choosing N and tq 
large enough. For the single bad component Qb of Q we cannot establish outright 
smallness. However, we will show that for a large enough n we have 

(6.33) WiTQbY^Wx^ xx^ -.x^ < 1 

Combining (|6.32p and (|6.33p it follows that for large enough N and tq we have 

ll(rQ)"llj>s:o"xJ>s:"^x,f xx" < 1 

This ensures the convergence of the series in ()6.30|) in X^ x X^ . Given the bounds 
in Lemmas 16.316.41 the proof of Proposition 16.21 is concluded. It remains to show 
that Q admits a decomposition which satisfies (|6.32p and (|6.33p . 
We begin with the easiest part, namely 

which will be included in Qg. Since the kernel K^e is bounded and rapidly decreasing 
at infinity we obtain 
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1 

which yields a r^c+i) gain, 

(6.34) W^JCcedrXoW 

N+ 2(13 + 1) ^ ll-^ollx™ 

For the second part uj^TZceXo of QceXo such a simple bound no longer suffices, 
and we need to use some cancellations. The final result is somewhat similar to the 
one above, in that it gains a power of r provided that /3 > 3/2. 



Lemma 6.5. The following estimate holds 



(6.35) WTcj^UceXoW 20-3 < \\xo\\xN 



Proof. Suppose that 

l|2;o||x„" = 1 

We denote g — uPlZceXQ and x ~ T^g. As before we also introduce the auxiliary 
variables y = p^x and h = p^g. The kernel Rce of TZce is bounded, rapidly 
decreasing at infinity and has symbol-like behavior at both and infinity. Then for 
the function h we directly estimate 

|t^-7^(1 +0/i(r)| + |r^+fT^(l + e)(9r - 2uj^d^)h(T)\ < \\xo\\x- = 1 
As in the proof of Lemma 16.41 we have 



yir, e(T)) = J U{t, (7, ar))h{a, C(fT)) da, 

where 

U{t, a,0 = sin (^C(r)H(r) \-\u)du^ 

Hence for y we use (I6.2ip amd (|6.22p to obtain the pointwise bound 

\yir, ar))\ < ar)-^ min{l, ^(r)^c^(r)-i} ^^-^+,^ (i + ^(^))-i da 
which we rewrite in the form 

(6.36) c.(r)|y(r,OI<r^min{l,^^c.(r)-i} / . da 
To bound we integrate by parts, 

y{r,i{r))^i{T)-^ h{a,aa))d, [\- cos (^i(T)"^ \{t) j\-\u)d,}j^ da 



^(^y (^cos (^^{T)iX{T) l\-\u)du^ -?jdMa,a'j)) 



da 



Estimating either |1 — costij < 1 or |1 — cosi;| < |f | this leads to a bound which is 
weaker than (|6.36[) . namely 

(l + gA2(r)A-^(a)) ^^ 

In a similar manner we evaluate D^-y, 

Dry{T,ar))^ J^^ ^^h{a,^{a)) COS (^e(r)U(r) ^\-i(w)du^ da 

= ^{r)-^ sin(^e(T)H(r)^ X-\u)du^ d.hia,^'^)) da 
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which leads to the same bound as in ()6.37p . Summing up, for z as in (I6.23P we 
obtain 



(6.38) W.,OI<r^min{U^.(.)-}y^ (l+e^^(.)A-^(.)) -^-- 

It remains to evaluate the integral on the right. If ^ < 2 then we can neglect the 
first factor in the denominator of the integrand and evaluate 

,.00 -AT-fci 

/ 1 !!! rf^<^-jv+^ 

X (l+a^(r)A-2(a)) ' 
However, if ^ > 2 then this factor yields rapid decay when 

i\\r)\-\a) > 1 

which corresponds to 

'^<T+(logO^+' 

Thus we obtain 



■da<(T + (loga''+i)-^+/^, ^>2 



(l-|-eA2(r)A-2(a)) 
Summing up, for z we have obtained the pointwise bound 

[ r^(T + (iogO''+')"'^+^ e>2 

This allows us to estimate norms, namely 

respectively 
Finally we obtain 

II^IIl- , L2n/~ L2 < 1 

N 2 JV 2 1/m 

20 + 1) 2(/3+l) 

and the conclusion of the lemma follows. 



□ 



Next, wc turn to the term Q^cX given by 

We will prove that QecX can also be included in Qg. The kernel RedO of TZec is 
bounded and decays rapidly at infinity. Then the contribution of the first term is 
easy to estimate using the L^L'^ type bounds on x and £_^x, 

(6.39) \\uJ^Recx\\ ^^^_s<\\x\\x«, S>0 

with S arbitrarily small. The bound for the second term in QecX is similar: 
Lemma 6.6. For 6 > we have 

(6.40) \\ToLjlCecDrX\\ „ + < \\x\\xN 



x„ 
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Proof. Set y = p^x. The solution xq = T^ujKecX is represented as 

/OO /"OO 

Uo{t, a)ujK,,{^) {d, - 2w{^di + 1)) y{a, d^da 
Integrating by parts we obtain 

POO poo 



/OO POO 
-dMr, aMa)K,,{Oy{a, ^)+2Uo{t, a)Lj\a)y{a, O^d^K^^iO d^da 

/OO poo 
^{cj){l+^)-^\x{a,£,)\d^dcj 



It Jo 

Hence 

*>CXD /'C30 



therefore 

A similar computation yields 



which leads to 
The desired conclusion follows. 



□ 



Finally we consider the expression QccX which has the form 

(6.41) QccX = oj^TZccX — 2ujJCoDt 
The first term is better behaved and can be included in Qg-. 
Lemma 6.7. For 5 > Q we have the following bound: 

(6.42) \\oJ^nccx\\^^^^^_, < \\x\\xN 

Proof. By the definition of the and spaces, it suffices to show that 

\\'T^ccx\\Ll<r\U^xUi+T-MA\Ll) 
This in turn follows by duality and dyadic summation from the bound 

(6.43) \\xw^)Kcf\\Ll < min(/l^-^ 1)|1/|U. 

For this we need to prove that the operator 7?.*^ is quasi-smoothing according to 
the following definition: 

Definition 6.8. A bounded operator T : Lp(M+) L^(R+) is quasi-smoothing if 
for each S > there exists Cs > so that 

(6-44) \\x[o,h) Tfhl < Cs mm{h-2-', l)\\fU2 

for all h>0. 
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We remark that the quasi-smoothing operators form an ideal under composition 
from the right within the algebra of bounded operators. Hence, given the expression 
of TZcc, it sufhces to show that the following operators are quasi-smoothing: 

^0, [^d^:l^o], K^ceK^ec, Jcc 

Recall that 

Jo i-v 

where 

\F{i, ,7)1 < min + 7y, + vYH^ + 1?^ " W""] 
LetFi(^,7?) :-p(0F(C,7?). Then 

For the first operator on the right-hand side one has 
(6.45) |^^i(e,r?)|X[i6[i,2]] <niin(e,l) 

which implies that the corresponding operator is quasi-smoothing, see the proof of 
boundedness of /Co in the previous section. For the second operator, observe 
that 

^3^X[i^[i,2]]<niin(l,(C + ^) ) 
by the rapid off-diagonal decay of F . Hence, 



sup 



X[i0[i,2]] ^i'^)dr\ < 11/11 J 



and therefore 

as desired. 

For the commutator [^d^,ICo] we have 

JO ^ ^ '? 

and one argues as before. 

The kernel of operator JCce^^ec is p{^)Ke{(,)Ke{ri), and is bounded by (1-f "(1+ 
rj)""". The quasi-smoothing property easily follows. Finally jTcc is quasi-smoothing 
due to the kernel bound (|6.7p . 

□ 

It remains to consider the second part of Qcc namely the expression LoJC^DrX. 
Since the kernel for /Cq decays at and at infinity, it is easy to establish the bound 

II^OIU^ < 1 

p/m pm 

It follows that 

(6.46) \\u;JCoDrx\\YN < \\x\\xn 

The difficulty is that there is no smallness in the above relation, and it is not 
possible to gain any smallness by letting r be large enough. To deal with this we 
reiterate: 
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Lemma 6.9. Suppose that n is large enough. Then 

(6.47) \\{TujlCoDrTx\\x» ^\\x\\x«- 

Proof. By (|6.46p and Lemma [6.41 it suffices to prove that for large enough n, 

(6.48) UcjICoDrTrgWi^Ll^ « \\g\\i^Ll^ 
We divide the operator JCq in two parts, 

with kernels 

<(C,r7) =ifo(C,r7)X[|i-i|<i], ^o'(C,'/) -^o(C,r/)X[|i_i|>i] 

The contribution of JC^'^ is non-resonant, we and we expect to gain powers of r 
from oscillations. Precisely, we will prove that 

(6.49) \\{uK^DrT){u;ri''DrT)g\\i^ <„ Mi^li^ 

Here the implicit constant depends on n, but that is not important since we gain a 
power of r. 

Assuming (|6.49p holds, in order to prove (|6.48p it remains to show that for large 
n we have 

(6.50) II {uKtDrTTgWi^Ll^ « bllz-L^^ 

Proof of (|6.50p : For another small parameter e to be chosen later we further 
divide /Cq into three parts, 

'^O ^ '^04 + ^^0,2 + '^0,3 

with kernels 

<'l V) - k<eK^{^, V), <'3 (C, V) = l?>e-<(^, V) 

The center part /Cq'2 enjoys better localization, while the two tails /Cq'i ^^'^ ^0 3 
are small. Precisely, 

(6.51) WujlC^llDrTgWiaa^^ + \\u;IC'^-'lDrTg\\i^L2^^ < e^bllz-L^^ 
It is easy to see that due to the supports of the kernels we have 

IC^ofiDrTJCt]f^-^ = 

and 

„(i,e(l+i) „ „„d,e(l+i) _ 
'^0,2 -^r J ^-0,3 = U 

Hence we obtain the decomposition 
(w/C^D^T)" = (tj/Co'jD^T)" 

n 
n 

+ 5](c^/c[Jfz?.r)^--i(c.<3^i?.r)(c./c^;^i?.T)"-^- 

l<j<k<n 
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For the middle part we will prove the bound 
(6.52) 

Combining this with (j6.5ip we obtain 



||(«CS:aD,T)'s||,-^;., < M!£M^||g||,.,,__ 



\\9\\l^L%„ 



Choosing e — n ^ this gives 



for a new constant C. Thus (|6.50p is established for n sufficiently large. 
We return to prove (|6.52p . Since 

d(Tidr]QUjJC'^'l{^, f]o)V{T, (Ti, ?7o) 



we can write the function 



in the form 



d(T2dr]iLjlCQ'2 [ Vo 



A2(ai) 



771 y(cri,cr2,?7i) ■ 



^Vn^i V{an-i,crn,rin-i) 



, ^d.e ( A^K-l) \ . A(t) 



A((T„ 



In order for the above integrand to be nonzero we must have 



f?/cA^(CTfc) 



- 1 



1 

< 

n 



1-1 
^70 



1 

< -. 

n 



rjk+iX^icTk+i) 
This implies that 

?7„A2(a„) <3^A2(t) 
Since £ < cr„, ^ < e^^ it follows that 

A2(a„) < Se^X'^ir) 
If T is sufficiently large this implies that 

Using the boundedness of /Cq 2 and of the transport along the flow (as \V\ < 
it follows that 



(r)|U. <m2(e)(Cc.(r)) 



/ dai j da2--- \\h{an)\\L2^^^^ da, 
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Changing the order of integration this yields 



\\y{T)U.<m\e)- 



(n- 1)! 



J (t - cr„ 



p/m 



Since 



cr(T) -i 

n 



we finally obtain 



„ „ 2. ^(C|loge|'^+^)'^,, „ 



Thus ([^3^ is proved. 

Proof of ([09l) : Denoting x = T{ujlCl^'^DrT)g, y — p^x and h = g we need 
to prove that 

(6.53) ||i?ry||i^ , < 

" 7JTT 

Due to the formula (|6.20p we have the integral representation 

poo \ f ^\ / ^ r(7 

■cosiri^X{a) / X~^{0)d6 ] y{a,r])dadr]ds 



where ^(s) — ^ ^^J^ and 77(5) = ?7^3^. In the support of the kernel ifp'^ we have 

> i. Thus the two oscillatory factors have 



1 



> - therefore 



ri2\(a) 



different frequencies, and we can gain if we integrate by parts with respect to s. 
Denoting 



u{s)^^^x{T)j x-'ie)de, v{s) ^ x{a) j x-'{e)de 

we write 

2 cosu(s) cosu(s) = cos(ii(s) + v{s)) + cos{u{s) — v{s)) 

We change the order of integration in the above expression for y and integrate 
by parts with respect to s. Since 

as 

we integrate the cosine and differentiate the rest to obtain 

rooroor^^^^^ 1 1 d 



mo \ / 



sin(u(s) + v{s))h{(7, rf) dsdrjda 



00 pOQ 



± 



A(r) 



It Jo 

noo poo 



A(r) 1 



A2(r) 



sin(D(r))/i((T, 77) dr/dcr 



r Jo 



,7) sin(7.(a))/i(a, 77) dr/da 
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We have 



Due to Theorem 15.11 the kernel Kq''' is bounded and decays rapidly at infinity 
therefore we can bound it by 

1 



We also have 



1 



(1 + 0(1 + ??) 

Hence the following rough bounds are valid: 



^2 ± 7^2 



^2 ± 7y2 



1 



4 2 7^2 



Inserting this in the bounds for Dry we obtain 



\DryiT,0\ <n 



oo poo per 



,2.^A(r) 1 



\h{a, r]) \ dsdrjda 



OC /'OC 



coir) 



A(r) 



T Jo A(cr) ^l,y(T)5 A2(t) 

A(r) 1 



|/i((T, r}) \ dr]da 



OC poo 



w(cr)- 



It Jo A(cr)Ocr)2r7 
This can be rewritten in the form 



^^\h{(7,r])\ d-qdu 



OC poo 



r^\Dry{T,S)\<n 

Taking weighted norms we obtain 

< -.2 



dsdrjda + 



rj2 

1 

TJ2 



oo /'OC 



t ^\Hovr])\ 
uj[T) J drjda 



drjda 



Then by Cauchy-Schwarz 



and further 



\DrV{T)h 



\Dry\\ 



2 ^1/ 



Thus (j6.5ip is proved, and the proof of the lemma is concluded. 
Proposition 16.21 follows. 

The proof of Proposition [6T] is also concluded. 



□ 
□ 
□ 

We now turn to the proof of Theorem 12.21 in Section [2] which estimates forward 
solutions e of the equation (12. 2|) . which we rewrite as 



Poe = f, Po = -df + dl + -dr + \{\ - 3Q(A(0r)2 - %Q(\(t)r)v^^) 
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under the action of the invariant vector field S — tdt + rdr ■ So far we have proved 
the bound p.3|l for e. In order to prove (|2.4p we write an equation for Se, namely 

PoSe^Sf+[Po,S]e 

A direct computation yields 

[Po, S] = 2Po-V, ^S{3Q{X{t)rf ~ 6Q{\{t)r)vio) 

Hence 

PoSe ={S + 2)/ + Ve 
A direct computation shows that 

' ' r2 (1 + i?2)2 ^ ^ 

Hence applying (|2.3p we obtain 



Then ([^^ follows since 
We remark that this requires 



iVo > iVi + 2/3 + 1 

The proof of (|2.5p is similar. 
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